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INTRODUCTION 

Considerable  effort  is  required  to  prepare  the  input  data  required  by 
the  Douglas  Three-Dimensional  Potential  Flow  Program.1  Extensive  hand 
calculations  are  necessary  to  provide  the  body  coordinates,  appendage,  body 
intersection  points,  and  appendage  coordinates.  These  calculations  are  then 
organized  into  the  Douglas  format  and  cards  are  punched  and  manually 
checked.  Because  this  procedure  is  so  time  consuming,  a computer  program 
package  was  developed  to  calculate  and  format  the  required  three-dimensional 
coordinates. 


The  development  of  this  package,  which  consists  of  three  computer 
programs,  is  complete.  The  programs  calculate  the  three-dimensional 
coordinates  for  axisymmetric  bodies  with  and  without  appendages,  calculate 
the  radial  and  angular  distribution  of  off-body  points  required  for 
propeller  thrust  deduction  calculations,  and  provide  punched  card  output  in 
the  format  required  by  the  Douglas  program.  These  programs  reduce 
coordinate  data  preparation  time  from  1 week  or  more  to  1 or  2 days,  making 
it  possible  to  provide  three-dimensional  potential  flow  analysis  in 
approximately  3 days. 


This  report  documents  the  procedures  for  using  the  computer  programs. 
The  text  is  divided  into  the  following  sections:  (1)  a description  of  the 
three  programs,  (2)  a definition  of  terms,  (3)  a specification  of  input 


quantities,  (4)  the  calculated  coordinates,  (5)  the  Douglas  program  format, 
C6)  off-body  points,  (7)  a description  of  the  schemes  used  to  input  the 


1 Douglas  Aircraft  Company,  Report  ES40622,  "Calculation  of  Nonlifting 
Potential  Flow  about  Arbitrary  Three-Dimensional  Bodies,"  by  J.  L.  Hess  and 
A.M.O.  Smith,  March  1962. 
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data,  (8)  typical  program  output  including  the  correspondences  between 
conventional  symbols  and  their  computer  representation,  and  (9)  a complete 
FORTRAN  listing  of  the  programs  and  the  nonsystem  subroutines  required  by 
them. 

Program  development  was  carried  out  on  a UNIVAC  1110  computer  system; 
however,  standard  FORTRAN  IV  was  used  so  that  the  programs  would  be  easily 
adaptable  to  other  large  computer  systems.  Although  the  FORTRAN  listings 
are  included  for  completeness,  a copy  of  the  programs  may  be  acquired  by 
contacting  the  author  (a  listing,  card  deck,  and  check  solutions  will  be 
provided). 


DESCRIPTION  OF  COMPUTER  PROGRAMS 
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The  three  computer  programs  described  in  this  report  calculate  and 
format  the  three-dimensional  coordinates  for  axisymmetric  bodies  with  and 
without  appendages.  The  programs  have  been  designed  to  require  a minimum  of 
input  data  and  to  provide  flexibility  in  the  specification  of  appendage 
shapes  and  locations  and  in  the  number  of  output  stations  along  the  body. 

Each  program  is  limited  to  a maximum  of  1000  output  points.  The  three 
computer  programs  are  as  follows. 

PROGRAM  1 - AXISYM 

] 

This  program  computes  the  three-dimensional  coordinates  for  an 
axisymnetric  body  without  appendages  and  with  one,  two,  or  three  planes  of 
symmetry.  The  final  results  consist  of  a maximum  of  100  equally  spaced  body 
stations. 

PROGRAM  2 - APNDG1 

E I 

The  second  program  computes  the  three-dimensional  coordinates  for  an 
axisymnetric  body  with  two  or  three  planes  of  symmetry  and  an  appendage 
configuration  located  in  the  afterbody  region.  Some  flexibility  is  provided 
for  the  number  of  output  stations  along  the  body  and  the  appendages. 

PROGRAM  3 - OFFBDY 

The  final  program  calculates  the  radial  and  angular  distribution  of 
the  off-body  points  required  for  the  construction  of  the  potential  flow 
velocity  profiles  needed  to  calculate  propeller  thrust  deduction.  The 
, starting  points  and  the  radial  and  angular  increments  are  program  inputs. 
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DEFINITION  OF  TERMS 


h 

h 

~C~ 
h ) 

1TJREF 


m,n 

t 

t 

~1T 

_C_;REF 


X-Y-Z 


AOB 

XS 

*0 


yz 

0S 


Appendage  chord 

Appendage  half  thickness 

Appendage  nondimens ional  half  thickness 

Value  of  -ji-  for  the  reference  appendage  half  thickness 
distribution 

Defines  a position  in  the  Douglas  formatted  coordinate  array 

Maximum  thickness  of  the  appendage  section 

Ratio  of  appendage  section  maximum  thickness  to  chord 

Value  of  for  reference  appendage  section  thickness 
distribution 

Cartesian  coordinate  system 
Appendage  nondimensional  chordwise  coordinate 
X location  of  the  off -body  points 
Location  of  the  appendage  leading  edge 
Appendage  leading  edge  offset 
Angular  increment  - off-body  points 
Radial  increment  - off-body  points 
Starting  angle  - off -body  points 


Subscripts 

i Point  identifying  the  beginning  of  the  nonredundant  portion 

of  a body  having  three  planes  of  symmetry 
j Arbitrary  point 

n Last  point 

o First  point 

R Value  at  the  appendage  root 

T Value  at  the  appendage  tip 
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SPECIFICATION  OF  INPUT  QUANTITIES 


BODY  CONTOUR 

A rectangular  cartesian  coordinate  system  (X,Y,Z)  is  used  to  describe 
both  the  input  and  output  coordinates.  Figure  1 shows  a typical 
axisymmetric  body  described  by  the  relationship  Z = f(x).  For  bodies  with 
one  or  two  planes  of  symmetry,  the  body  contour  is  described  by  an  array  of 
points  from  (XQ,  ZQ)  to  (Xn,  Zn)  (figure  1A).  For  three  planes  of 
symmetry,  only  the  nonredundant  portion  of  the  body  is  input;  therefore,  the 
body  is  described  by  points  from  (X^,  Z^)  to  (X^,  Z^)  (figure  IB).  Regions 
of  large  curvature  should  contain  a sufficient  nunber  of  points  so  that  the 
contour  will  be  accurately  represented. 

APPENDAGE  CONFIGURATION 

The  appendage  thickness  distribution  at  all  radial  stations  is  derived 
from  the  reference  nondimens ional  half- thickness  distribution 
= f(XD)  input  to  program  2.  The  chordwise  variation  of  the  nondimensional 

D 

half  thickness  at  the  appendage  root 

_h_)  t/C)R  f(x  ^ 
c R t'/OpEF  B 

and  at  the  appendage  tip 

= t/C)T  f(X  ), 

t7Cr ref  b 

is  used  to  calculate  a linear  variation  from  the  tip  to  the  root  of  the 
appendage  thickness  t.  Figure  2 shows  a typical  appendage  cross  section; 
however,  appendages  are  not  restricted  to  airfoil  shapes.  The  additional 
input  appendage  quantities  are  defined  in  figure  3.  Also,  to  provide 
flexibility  in  program  2,  any  appendage  shapes  found  in  figure  4 may  be  used. 
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Figure  2.  Typical  appendage  cross  section. 


CALCULATED  COORDINATES 


APPENDAGE  COORDINATES 

The  calculated  appendage  coordinates  are  described  by  section  numbers 
corresponding  to  various  Z locations  (figure  5).  Each  section  contains  the 
(X,  Y)  coordinates  describing  the  appendage  cross  section.  Only  that 
portion  of  the  appendage  parallel  to  the  positive  Z axis  is  calculated, 
while  the  coordinates  parallel  to  the  positive  Y axis  are  obtained  from 
symnetry. 

BOD'.'  COORDINATES 


Any  body  section  not  containing  an  appendage  is  represented  by  10 
equally  spaced  points.  For  one  and  two  planes  of  symmetry  the  body  section 
is  divided  into  20-degree  increments,  and  for  three  planes  of  symmetry  the 
body  section  is  divided  into  10-degree  increments.  For  body  sections 
containing  an  appendage,  the  10  body  points  are  distributed  equally  with  the 
first  and  tenth  points  coinciding  with  the  appendage  intersection  points 
(figure  5). 


(X.Y)  APPENDAGE 
POINTS  


Figure  5.  Calculated  coordinates  - body  and  appendage. 


DOUGLAS  THREE-DIMENSIONAL  COORDINATE  FORMAT 


The  Douglas  Three-Dimensional  Potential  Flow  Program  requires  the 
input  body  and  appendage  coordinates  to  be  arranged  in  an  array  of  (n,  m) 
points,  where  n identifies  the  column  and  m identifies  the  position  in  the 
column.  The  n and  m designations  follow  the  convention  shown  in  figure  6. 
An  observer  located  in  the  flow,  oriented  so  that  the  m values  increase 
upward,  also  sees  n values  increasing  to  the  right.  The  Douglas- formatted 
output  is  arranged  so  that  the  body  sections  are  followed  by  the  appendage 
sections  located  in  the  positive  Y-Z  quadrant. 

The  Douglas  input  scheme  also  requires  the  use  of  a status  word  to 
indicate  the  beginning  of  new  sections  and  to  identify  the  last  data  point. 
The  input  data  cards  must  also  contain  sequence  nunbers  arranged  in 
ascending  order.  The  cards  punched  by  programs  1,  2,  and  3 have  the 
following  format: 


CARD 

COLUMNS 

1-10 

11-20 

21-30 

31 

32-41 

42-51 

52-61 

62 

77-80 

VARIABLE 

X 

Y 

Z 

STAT 

X 

Y 

Z 

STAT 

SEQ 

FORMAT 

F10.5 

F10.5 

F10.5 

11 

F10.5 

F10.5 

F10.5 

11 

14 

SEQ  is  the  sequence  number,  and  STAT  (the  status  word)  is  defined  as  follows: 

STAT  = 1 - data  point  is  the  beginning  of  a new  n line 

STAT  = 2 - data  point  is  the  beginning  of  a new  section 

STAT  = 3 - data  point  is  the  last  point. 

The  status  word  is  left  blank  otherwise. 


z 


n COLUMNS 


t 


Figure  6.  Organization  of  input  points  into  rows  and  colunns. 
(Adapted  from  Douglas  Aircraft  Company,  Report  ES40622, 
"Calculation  of  Nonlifting  Potential  Flow  About 
Arbitrary  Three-Dimensional  Bodies," 
by  J.  L.  Hess  and  A.  M.  0.  Smith, 

March  1962.) 


The  method  of  inputting  coordinates  to  the  Douglas  program  also 
requires  a change  in  the  coordinate  system  origin  for  bodies  with  three 
planes  of  synmetry.  In  this  case,  the  coordinate  system  origin  is 
translated  to  (X^,  (figure  IB)  so  that  the  nonredundant  portion  of 
the  body  begins  at  X = 0.  The  calculated  body  and  appendage  coordinates  are 
adjusted  by  the  relationship  Xj)^^  = X^q^  - X^.  The  origin  change  is 
performed  within  programs  1 and  2 and  is  reflected  in  the  Douglas- formatted 
output.  The  coordinate  system  origin  remains  the  same  for  bodies  with  one 
or  two  planes  of  synmetry. 
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OFF- BODY  POINTS 


Program  3 computes  the  radial  and  angular  distribution  of  off-body 
points  required  in  the  construction  of  potential  flow  velocity  profiles. 
The  starting  points  X^g  and  Zg  and  the  radial A^z  and  angular  aq 
increments  are  defined  in  figure  7.  Hie  punched  card  output  of  program  3 
has  the  same  format  as  described  in  the  previous  section.  It  should  be 
noted  that  the  coordinate  system  change  in  origin  for  bodies  with  three 
planes  of  synmetry  must  also  be  applied  by  the  user  to  the  starting  XQB 
values  of  the  off -body  points. 


INPUT  TO  COMPUTER  PROGRAMS 


METHOD  OF  INPUTTING  TABLES 

The  input  tables  of  body  and  appendage  coordinates  are  described  by 
functional  relationships  of  the  type  Y = f(x).  All  tables  of  this  type  are 
input  according  to  either  scheme  1 or  scheme  2. 

SCHEME  1 


F Format 


CARD 

COLUMNS 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

X. 

l 

Yi 

X.  , 
l+l 

Yi+1 

X.  , 
1+2 

Y.  , 

i+2 

SCHEME  2 


F Format 


CARD 

COLUMNS 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

X1 

X2 

• • • 

X3 

"n 

• • • 

. . . 

Y1 

Y2 

• • • 

Y3 

Yn 

• • • 

• • • 
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DEFINITION  OF  FLAG  CARD  VARIABLES 

The  flag  card  variables  conmon  to  programs  1 and  2 are  described  below. 
All  variables  are  of  the  I format  and  must  be  right  justified. 

Variable  I Description 


If  IREAD  = 1,  all  data  tables  are  input  according  to 
scheme  1 

If  IREAD  = 2,  all  data  tables  are  input  according  to 
scheme  2 

If  IFRINT  = 0,  only  the  input  data  and  the  Douglas 
formatted  output  are  printed 

If  IFRINT  = 1,  the  input  data,  intermediate  results, 
and  Douglas  formatted  output  are  printed 

If  IPLOT  = 0,  no  plotting  is  performed 

If  IPLOT  = 1,  a three  dimensional  plot  of  the  body  is 

made 

If  I PUNCH  = 0,  no  cards  are  punched 

If  IFUNCH  = 1,  the  Douglas  formatted  ouput  is  punched 

If  IDOUG  = 0,  the  Douglas  formatted  output  will  not  be 
calculated 

If  IDOUG  = 1,  the  Douglas  formatted  output  will  be 
calculated 


t 


The  number  of  planes  of  synmetry 


The  following  defines  the  additional  flag  card  variables  used  in 
programs  1,  2,  and  3. 


PROGRAM  1 

Variable 


Description 


INPUT  SCHEME  FOR  PROGRAM  1 


LABEL  CARD 


Card 

Input 

Columns 

Format 

Quantity 

1-80 

A 

Label 

Alphanumeric  characters 

FLAG  CARD 


Card 

Columns 

Format 

Input 

Quantity 

1-3 

I 

IREAD 

4-6 

I 

IPRINT 

7-9 

I 

I PLOT 

10-12 

I 

I PUNCH 

13- IF 

I 

NSW 

NSYM  = 1,  2,  or  3 

16-18 

I 

I DOUG 

19-21 

I 

NPTS 

20  4 NPTS  4 200 

22-24 

I 

NB0DY 

NBODY  4 100 

TABLE  1 - Z = f(X) 


Maximum  number 
Minimum  number 


of  entries  this  table  may  contain 
of  entries  this  table  may  contain 


is  200. 
is  20. 
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INPUT  SCHEME  FOR  PROGRAM  2 


LABEL  CARD 

Card 
Col  unns 


Format 


Input 

Quantity 


1-80 

A 

LABEL 

Alphanumeric  characters 

FLAG  CARD 

Card 

Input 

Colwms 

Format 

Quantity 

I PRINT 


I PUNCH 


NSYM  = 2 or  3 

20  < NBODY  < 200 
5 < NFIN  < 25 

If  I DOUG  = 0,  IFWD  <60 
If  I DOUG  = 1, 
and  NSYM  = 2 , IFWD  < 40 
NSYM  = 3,  IFWD  < 30 
If  I DOUG  = 0,  I AFT  < 25 
If  I DOUG  = 1, 
and  NSYM  = 2,  I AFT  < 10 
NSYM  = 3,  I AFT  < 20 
5 < IFIN  < 14 


TABLE  1 - Z = f(X) 


Maximum  number  of  entries  this  table  may  contain  is  200. 
Minimim  number  of  entries  this  table  may  contain  is  20. 


CONSTANT  CARD  1 


Card 

Columns 

Format 

Input 

Quantity 

1-10 

F 

CT 

11-20 

F 

°R 

21-30 

F 

xs 

31-40 

F 

X 

41-50 

F 

o 

ZT 

51-60 

F 

ttjref 

61-70 

F 

t ) 

-rr}  t 

71-80 

F 

~tjr 

TABLE  2 - -£-)REF  = f(X)„ 


Maximun  number  of  entries  this  table  may  contain  is  25. 
Minimun  number  of  entries  this  table  may  contain  is  5. 


24 


INPUT  SCHEME  FOR  PROGRAM  3 


LABEL  CARD 


Card 

Input 

Colunns 

Format 

Quantity 

1-80 

A 

LABEL 

Alphanumeric  characters 

CONSTANT  CARD 


Card 

Columns 

Format 

Input 

Quantity 

1-3 

I 

NPTSZ 

See  note  below 

5-7 

I 

NTHETA 

11-20 

F 

^B 

21-30 

F 

ZS 

31-40 

F 

A 

yz 

41-50 

F 

es 

51-60 

F 

4e 

In  degrees 

NOTE:  Repeat  the  constant  card  for  each  desired  value  of  Xgg  (.the  total 

number  of  points  <1000.).  To  initiate  the  output  of  the  complete  set 
of  calculated  body  points,  a constant  card  with  NPTSZ=0  must  be  used. 


COMPUTER  PRINTOUT 


INCLUDED  ITIMS 

All  input  quantities  to  programs  1,  2,  and  3 are  printed  out  to  allow 
rapid  checking  for  errors.  The  results  of  the  computer  solutions  are 
tabulated  and  identified  with  the  computer  correspondences  to  conventional 
symbols  given  in  this  section.  The  following  is  a list  of  the  items 
included  in  the  output  of  each  program. 

PROGRAM  1 

1.  Control  parameters 

2.  Input  body  coordinates 

3.  Calculated  three-dimensional  coordinates 

4.  Calculated  coordinates  - Douglas  format 

PROGRAM  2 

1 . Control  parameters 

2.  Input  body  coordinates 

3.  Input  appendage  data 

4.  Intersection  point  summary  - appendage  and  body 

5.  Calculated  appendage  coordinates 

6.  Calculated  body  coordinates 

7.  Calculated  coordinates  - Douglas  format 

PROGRAM  3 

1.  Control  parameters  for  each  X location 

2.  Complete  set  of  off-body  points  in  the  Douglas  format 
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PRINTOUT  CORRESPONDENCES 


SYMBOL  COMPUTER  REPRESENTATIVE 


REF 


B 


Y 

Z 


h. 

zs 

ZT 

Ay, 

Ae 

es 


TYPICAL  COMPUT®  RUNS 


CORDR 

CORDT 

H 

HC 

TCREF 

TCROOT 

TCTIP 

X 

XB 

XOFSET 

XSTART 

Y 

Z 

ZROOT 

ZSTART 

ZTIP 

DELYZ 

DELT 

TSTART 


Several  examples  have  been  included  to  illustrate  the  results  produced 
by  programs  1,  2,  and  3.  Figure  8 shows  the  results  of  program  1 for  an 
axisymnetric  body  with  two  planes  of  synmetry,  while  figures  9 and  10  demon- 
strate the  results  of  program  2 for  the  cases  of  two  and  three  planes  of 
symnetry.  Hie  computer  printouts  for  the  three-dimensional  coordinates  and 
off-body  points  for  the  example  shown  in  figure  10  are  found  in  appendix  A. 


APPENDIX  A.  COMPUTER  PRINTOUT  EXAMPLES 


calculation  the  three  dimensional  coordinates 

FOR  »N  AxISymMETrIC  BODY  WITH  APPENDAGES 

having  j planes  or  symmetry 


TEST  CASE  R 


CONTROL  PARAMETERS 


IREAD  > 

2 

IpRlNT  ■ 

1 

I PLOT  ■ 

1 

ipunch 

NBODY  - 

SR 

nfin  * 

17 

IFWD  ■ 

30 

IaFT 

nsym  . 

3 

(DOUG  ■ 

I 

1 F 1 N . 

12 

mis  PAGE  IS  BEST  QUAIlT*  PRACTICABLE 
KBOa  COPY  FUKUISH®  10  DOC  ^ 

INPUT  I'on''  COORDINATES 


X z 


\ 

. 750000o»07 

.*375000*01 

7 

• 300oOOo*OZ 

.*375000.01 

1 

. jsoonoo*o? 

.*375000*01 

M 

• 9000POo»0? 

.*375000*01 

«i 

. 9500000*07 

•A375000.pl 

A 

.500000o»o? 

^*375000.01 

7 

•550000o*07 

. *375000*0 1 

R 

. *000000*0? 

.*375000*01 

o 

. A50000o*0? 

.A375000.pl 

l n 

. 700000o»0? 

.*375000*01 

i 1 

. 750o00o*o? 

.*375000*01 

1 7 

. *00000o»07 

i A 375000*0 1 

1 I 

.*50000o*07 

.A375000.0I 

1 M 

. 900000o»0? 

. A375000*0 1 

1 S 

. 950000o*0? 

. *375000*0 1 

1 A 

• 1 000000*03 

.*375000*01 

t 7 

. 1 0S0O0q»03 

.*375000*01 

t A 

. 1 100000.03 

.A37SOOO.pl 

1 9 

. 1 |5onOo*OS 

.*375000*01 

?n 

. 1 1 9q00o*O3 

• A 3 7 5000*0  * 

71 

. 1 707000.03 

.*375000*01 

?7 

. 1 77| ?9p*oi 

.*370000*01 

7\ 

. 1 733*130.03 

.*33*000*01 

?*♦ 

. 1 7**5R7o.n3 

. *7*9000*0 1 

?«; 

. l?*771o.03 

.*730000*01 

?6 

.1?*9**0»03 

.*1*5000*01 

77 

. 1 7*7000.03 

.*093000*01 

?« 

.17991*0.03 

.*00*000*0' 

7 9 

. 1 3O*79o.01 

.5*91000*01 

30 

. 1 310*30.03 

.575*000*01 

3 1 

. 1 33o*7o»03 

.5597000*01 

3? 

. 1 3*77| 0.03 

.5*1 |000*0l 

33 

. 1 3*90*0.03 

.5707000*01 

3M 

. 1 3*7000.03 

.*95*000*01 

IS 

.13791 *0.01 

.**79000*01 

3A 

. 1 3*97*0.03 

.***9000*01 

37 

. 1 3973*0*01 

. * 1 9*000*01 

1 A 

, 1 909*0o»0 1 

• 3*73000*0  1 

39 

. 1 * 19*90.03 

.3*75500*01 

*♦0 

. 1 *7*B9o.01 

.317 1570*0 1 

. 1 97*390.n3 

. 79 | 8 7 *0*0 1 

97 

. 1 **?*9q.03 

.7*98730*01 

9 3 

. I**«790t01 

• 75 1 300O»O 1 

9«| 

. 1**3390.03 

.73**1 30*01 

9*; 

. 1 9*9390.03 

.7175300*01 

96 

.1**3390.03 

.7004*00*01 

97 

. 1 9A«3?o»03 

.1*375*0*01 

9 A 

. 1 977**0.01 

, 1 *990*0.01 

U«> 

. 1 *7***0. ni 

. 1 559000*01 

Sn 

, | 9*0*90.0  3 

, 1 9 | 5970*0 1 

*=>1 

. l9P599o.ni 

.17***00.01 

S2 

, | *9oA90.03 

< 10*0000*01 

S3 

. 1 *9SA9O.03 

. ’0*8*00*00 

sm 

. 1*00**0*03 

.7*00000*00 

XBIS  PAGEE  IS  BEST  QUALITY  PRACTICABLE 
JROU  OQPY  riMilSHED  10  DD,Q - 


INPUT  APPENDAGE  oat* 


start. 

, i ji?7io*oj 

xoeset* 

,4000000*00 

ZTIP* 

1 000000*02 

CORDT. 

, 400000°*0| 

CONOR* 

;RA7AOOO*OI 

Tc«fr. 

, 1 I 98g40*0o 

TCT IP. 

1 000000*00 

TCROOT. 

. 1 250000*00 

XB 

MC 

1 

,000000° 

.0000000 

2 

,■*94  1 00°-02 

. |57O900"0 1 

J 

. 1 Z50soO.O| 

, |89| JOO-OI 

4 

. 250 1 00°*° 1 

,2414700-01 

5 

.5002 | o°-0 | 

. 3555200-01 

A 

, 7t9?BoO-0| 

.4204300-01 

T 

. 999  3goO-0 l 

,4A8lI°0-0l 

A 

, | 499590*00 

.5343100-01 

9 

• 1 948  740*00 

.5735800-01 

10 

,7448970*00 

.5932200-01 

1 1 

, 7948 | 40*00 

.5944200-01 

1 2 

• 1498550*00 

, ^ 747400-0 1 

1 J 

.4947930*00 

.524 (400-01 

It 

.5498 J50*00 

.4557700-01 

1 5 

,*9477j0»00 

.3458500-01 

1 4 

• 7 997  | | 0«0p 

.2425100-01 

17 

• I00000°*0| 

.oooonoo 

EHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
ratJM  COPY  imtilSHSU  TO  DJDC 


lTE«SrCT|0N  POINT 

SUMMARY  « APPenDaGE 

AND  BODY 

» 

r 

2 

1 

i 1 3*127  |0*03 

, 0000000 

.591 | 153*01 

2 

. 1 3m7**TA.03 

. 3409397*00 

.5320249*01 

3 

. 1 357253.03 

»977?7**0*00 

.S27793g*n| 

M 

. 1 3*J  *93*03 

,5404454*00 

• 50  3 I A00*0 1 

s 

.t37i9*3*os 

,*297829*00 

.9*1*985*01 

* 

. ■ 3d i 1 *7*03 

.41 17277*00 

^45*81*5*01 

7 

* 1 39o’95.03 

.5*54299*00 

^9338751*01 

« 

. |9oo7*2»OS 

.9991913*00 

.90**81 i *01 

9 

. |9|O705*03 

. 9028007  *00 

.3742079*01 

10 

.1920779*03 

.2934429*00 

.3935219*01 

1 1 

. 1930*93*03 

.1412121*00 

.3|0|08**0l' 

1 2 

. M**  | 0*9.03 

.7939229-07 

,27*0940*01 
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IBIS  PACE  IS  BEST  QUALITY  PRACTICABLE 

Y80*  COPY  FURNISHED  TO  D<DC 


Anrr  N n A GE  SECTION 

X 

.13967(0*0 

MImBFR  | 2* 

y 

. i ooonoo‘i’2 

apPenoagF  SEC  1 1 On 

x 

number  ? z* 

Y 

.9926375*01 

. pnoonoo 

.139*210*03 

.0000000 

. 1 35  1 2 1 o*o3 

. JI66RS76*  JO 

. I 350757*03 

.2783236*00 

. 1 3*i  *•7  1 0*03 

•3515956*00 

. 1 355799*03 

. 3666572.00 

. 1 3 6 9 7 1 0 * 0 3 

.9306591*00 

. 1 369379*03 

.9991663*00 

.137371 o»o' 

.9500007*00 

. 1 373963*03 

.9693995*00 

. 1 3827 10*03 

.9352239*00 

. 1 302598*03 

.9539320*00 

.I3«17 | n»o3 

.3971910*00 

. 1 391*32*03 

.9192126.00 

. 1 9 oO  7 | 0*03 

. 39705*9*00 

. 1 9007 1 7*o3 

.3567601 *00 

■ion 9710*03 

•2799993*00 

. 1 90990  1 * 0 3 

.2862085*00 

. 1 9 | 8 7 | 0 ♦ o 3 

.(9*831 3*00 

, | 9 1 8996*03 

. 2052929.00 

.19277 | 0*03 

•106271 |*00 

. 1 927970*03 

. 1 108393.00 

. M.iA7io*03 

.5193707-07 

. 1 937q55*03 

.5369815-07 

APf-EMOAGF  SECTION 

NUMBER  3 2* 

,8852750*01 

appemoage  section 

v 

NUMBER  9 2* 

y 

.8779125*01 

• 1 3957 | 0*03 

. qoooooo 

A 

. 1 395210*03 

.0000000 

. 1 350799*03 

.2897995*00 

. | 399837*03 

. 301 2655*00 

. | 359979*03 

. 38 1 76«8*00 

, 1359969*03 

.3968803*00 

, | 3<,909B*o3 

. 9676789*00 

. I 363717*03 

.9861906*00 

. 1 37371 7*03 

, m“B6983»00 

. 1377971 *03 

.508032 1*00 

.1382386*03 

.9776905*00 

. 1 387229*03 

,9913991 *00 

.1391 555*03 

.93 1289 | *00 

.1391 977*03 

.9983557*00 

. 1 9(j07  2 9 *o3 

• 37  1 96  38*00 

. 1 90073  1 ♦ 0 3 

. 386 1679*00 

. 1 909893*03 

•2980877*00 

. 1 909989«q3 

. 3098869.00 

.191 9n*2*o3 

. 2 1 37539*00 

, 1 9 19238*03 

.2222199.00 

.1978231*03 

• 1 I5  9o75*oo 

.192099  1 ♦ 0 3 

. | 199757*00 

. 1 9 J 7 9QU*fj  3 

• 5S«S973-'J7 

. 1 937795*03 

.5807031-07 

A P P ( 1 1 0 A 5 E SECTION 

r H'.HER  5 7* 

. 7705500*01 

APPENEACif  SECT  ION 

number  6 z* 

.713 1875*0 1 

X 

. 1 399710*03 

. 0000000 

,1399710*03 

.0000000 

. 1 399379*03 

.3127169*00 

.1398921 *03 

, 3292079.00 

. 1 .35909.1*03 

.91  1 9?  | 9*(J.) 

. | 353633*03 

.9771035*00 

.136 3 386*07 

. 509 7 027  *U0 

. 1 363055. 03 

.5732199.00 

. 1 3 7 2 T 79  »|J3 

. 5773759., jn 

. 1 377977*03 

.5967 1 96.00 

. 1 397  "67*  '3  3 

. 5 1 015  7 7 .J'l 

.1381 900*03 

.5287662.00 

. 1 39  | 'ino*o3 

. 965 97 7 7 *UN 

.1391372*03 

. 9829988.00 

. 1 9 1 1 7'  7 3 1;  ♦ i.|  3 

.901871  1 « u 1 

, 1 9 fin  795*03 

. 9 1 55798.00 

. 1 1 1 "11  76  • 'J  3 

. 3 2 1 6 8 * 2 • J " 

. 1 9 1 1 ' 1 6 » « 0 3 

. 3339859.00 

. | 9 1 9 9 | 9*l|  .3 

. 2J06  75  9 *U0 

, 1 9 | 959Q.Q3 

.2391 369.00 

. | 928752*0 3 

. | 295939*00 

. 1 92901 2*03 

.1791121 *00 

. 1939090*03 

•60781 39-07 

.1938935*03 

.6299297-07 
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IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
mM  COPY  PUK&LSHED  TO  LDC  


r ■■  i p a ti r Section 

NUMBER  7 z« 

.6558250*01 

APPfNCA&E  SECTION 

number  0 z« 

13437 | 0*o3 

. 0°00000 

. 1347210*03 

, OOOOoOO 

1348963*03 

. 3 356784*00 

. 1 348006*03 

.3471493*00 

1 353? 1 7 ♦ 0 3 

.4472151*00 

. 1 352802*03 

, 95732*7*00 

1 36. '774*03 

• 54  1 7 7 7 0 ♦ 00 

. 1 362393*03 

.5602392*00 

1 3 7 27  3 1 ♦ 0 3 

•5660634*00 

.137 1985*03 

.5854072*00 

l 3^ l 7 

,5474748*00 

. 1 381576*03 

.5661833*00 

1391245*03 

.4995703*00 

,1391168*03 

.516641 9*00 

14  ij  0757*0  3 

.4302785*00 

.1400759*03 

,9999821 *00 

1 4 | 0259*03 

.3452846*00 

.141 0350*03 

. 3570838*00 

1 4 | 9 7 66*0  3 

.2475974*00 

. 1 4 1 9942*03 

. 2560584*00 

1 429773*03 

. 1 3 36802*00 

. 1 429533*03 

. 1 382484.00 

|M387P0*03 

.6470354-07 

. 1 4 39 1 25*03 

.6691462-07 

_'.na<«E  SECTION 

m ijm  p r si  v z = 

.5411 on0»0 1 

appenoa&e  section 

NUM0FP  10  2» 

X 

V 

X 

r 

13427|Q*o3 

. 0 0 0 0 0 0 0 

. 1342710*03 

• 0000000 

1 3 4 75  48*03 

. 3586203*00 

. 1 347524*03 

. 3592251  *00 

1382386*03 

. 4774387*00 

. 1 352342*03 

.4740502*00 

1 362062*03 

•578751 3*00 

.1361 »89*o3 

.5828 327*00 

1371738*03 

•6047510*00 

. 1371653*03 

.61 14281*00 

1381414*03 

•5848919*00 

.1381 337*u3 

.5938372.00 

1341090*03 

•53371 34*00 

.1391 042*03 

.5443504.00 

l 9q076**o3 

. | 40077 1*03 

.471 1710*00 

1410442*03 

•3688831*00 

.1410530*03 

. 3001  809*00 

14201 | B«03 

•76451 95*00 

. 1 420320*03 

.2 742338.00 

1 4 ? 9 7 9 4 ♦ 0 3 

. | 428  | 66*0(7 

.14301 44*03 

, | 989989*00 

1474473*03 

.6412571-07 

. 1 440001 *03 

.7253123-07 

: I !)  A <"■ F SECTION 

VlMlU  R )|  / = 

.5411 Q0O*D 1 

APff.  rjrAGl  SECTION 

number  i?  2* 

i 3 ♦ ? 7 i n ♦ o ' 

• n 0 0 p n P n 

. 1342710*03 

.noooooo 

1 v»  7 r>  n rj  ♦ j } 

• 359p2?9*(jo 

.1347476*03 

. 3604347.00 

1 15  2 2 9 7*03 

. 4 756671*0  <1 

. 1 757253*03 

,9772790*00 

1 7 6 1 4 | 6 ♦ fl  3 

. 5 8 6 4 1 4 0 . 0 ,7 

.1361843*03 

.59Q9?5M*0n 

1 3 7 1 5 h « ♦ o 3 

. 6 1 0 1 057*00 

.1371 483*03 

. 6 7 9 7 6 2 9 ♦ 0 0 

1 7 » 1 5 9 ♦ 0 3 

• A/’  ? 7 >'  ? 44  ♦ jri 

.13811 82*03 

.611 7777.00 

I 3 4 0 4 •>  3 ♦ 0 3 

, c,5  M 9 n 7 M *. ji 

. | 390945*03 

.5656244.00 

I w p • ! 7 7 7 ♦ ;l  1 

,4876661  *,•>'; 

. |Mnn7P7*o3 

.999191  3*00 

i ♦ | ’ M 7 *0  * 

, 1 404.:.  j,' 

.1410705*03 

. 9020007*0" 

1 470872*03 

. 283948 | *00 

. 1420724*03 

. 2936624.00 

1430493*03 

• 1 55080  3*00 

. 1430843*03 

.1612121 *00 

1 4 405  3 7*0 3 

.7593676-07 

. 1441064*03 

, 7934229-07 

.5984625*0 
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.5*Mi  noo*o  i 


.591  | 15  3*01 


c»lcul»tep  body  coordinates 


'■’s 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  COPY  FURNISHED  TO  DDC  


ROOT  SECTION  number  | X.  .2500000*02 


y 


7 


.0000000 
. I I07006»0| 
,?l*0377*0| 
, 3 1 8 7 9 9 * * 0 | 
, 1 0 * 7 7 A R * 0 | 
.9883530*0 | 
,5S?0’0R*0| 
.599o538*0| 
. 4 2 7 fl  | MR  * 0 l 
,4375000*01 


.4375000*01 
.4276150*01 
.5*9059 I *0  I 
• 5S709 | 3*0  I 
.1663534*01 
.1097775*01 
.3167505*01 
.2160365*01 
. I I 070  I 5 *0  I 
.6716854-05 


BOOT  SECTION  NUMBER  3 X.  .3337113*02 


boot  section  number  2 x-  .29ie794*02 


Y 

,0000000 

J I 107006*01 

,21 *0377*01 

.3167197*01 

.1097766*01 

jl863S3o*OI 

.5670909*01 

^6990536*01 

.6778116*01 

.4375000*01 


6375000*01 

6278119.01 

5990511.01 

5520913.01 

1663536.01 
1097775*01 

3187505.01 
2180385*01 

1107015.01 
8715855-06 


I 

] 

BOOT  SECTION  NUMBER  1 X«  .3766119*02 


Y 


2 


.oonoooo 

. I I 07004*0 | 
.21 8p377*0| 
,3167197*01 
.1097768*01 
.1663530*01 
.5570909*01 
•5990538»0 | 
,627b|18*0| 
.6375000*01 


.4375000*0) 
. 4778 | 19*0 i 
.5990511*01’ 
.55709 i 3*0  I 
.1883536*01 
,1097775*01 
.3187505*01 
.2180385*01 
. I 107015*01 
,6715855-05 


Jooooooo 

. I 107006*01 
,7)80377.01 
.31*7197*01 
| 1097746*01 
,98*3530*01 
.5570909*01 
,6990538*01 
,477* | 98*0  I 
^4375006*01 


,6375000.01 
,6278199.01 
,5990511 *01 
.5520913*01 
J 1683536*0 | 
,1097775*01 
,3187505*01 
J 2 I 60365*01 
, I 107015.01 
.8795856-05 


BOOT  SECTION  NUMBEP  5 X.  .9179825*07 


BODY  SECTION  NUMBER 


6 X.  .9593531*02 


Y 


2 


.0000000 
. I 107006*0  I 
,2I6o377»Oi 
. 3 I 67997*0 | 
. 9097748*0 | 
.9883530*01 
.55?p909*0| 
,S99o5  38*0 1 
,4?76|98«0| 
.4376000*01 


.6375000*01 
.427*1 99*01 
.5990591 *0 i 
.5520913*01 
.9863534*01 
.9097775*01 
. 3 I 67505*0  I 
.2180385*01 
. | 107015*01 
.8716855-05 


iooooooo 

. I 107004*0  1 
.21*0377.01 
j3l*7997«oi 
.9097778*01 
;i6*353o*OI 
t 56  70909*0  I 
j59905 38*01 
,4778196*01 
l 4376000*01 


4375000*01 

4778150*01 

5990591*01 

5520913*01 

9683536.01 

9097775.01 
3187505*01 

2180385.01 
1 107015*01 
8795855-05 


«on Y SECTION  NUmbER 

7 X. 

.5012236*07 

BOO Y SECTION  NUMBER 

8 X.  ,5930999*02 

y 

2 

Y 

2 

,0000000 
. I 107004*01 
,21*0377*01 
. II *7997*0| 

, 909,74*. 0 | 

. 9*6 15 10*0 1 
, 55  ?0909»0 1 
.5990538*01 
.677*196.01 
.4375000*01 


.4375000*01 
.6278 1 99»o 1 
.5990591 »0i 
.55709 i 3*01 
•96*3534.0) 
.9097775*01 
.3187505*01 
.7180365*01 
* I 107015*01 
.8795655-05 


Jooooooo 

. I I 07004*01 
J 2 I *037  7»o I 
.31*7997*01 
Ji0977A8*0l 
,96*3530*01 
,5570’09*0I 
,5990536*01 
J*?78|98*01 
• 4 375000*0  I 


4375000*01 
6778199.01 
5990591*01 
5520913*01 
9**3536*01 
9097775*01 
3187505*01 
2180385*01 
1 107015*01 
8795*55-05 


39 


«IS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
fROtf  COPY  FURBISHED  TO  DOC 


floor  SECTION  NUmbE 8 

9 X.  ,5899*50*02 

BOOT  SECTION  NUMBER 

M 

10  X*  .6768356*07 

* 

,0000000 

2 

•6375000*01 

T 

jooooooo 

* 

>375001*01 

. | 1 0 7006  *0 | 

.6278|89*q| 

, 1 107006*01 

.6278150*01 

, 2 1 8o37  7 *0 | 

.5990591*01 

>1  *0377*01 

j 5990592  *0  j 

,3I*7997»0| 

.5570913*01 

,3187998*01 

,55209l«*0l 

, 90977***0 | 

.9883536*01 

,9097769*01 

>883536*01 

, 9**3530*0 | 

.9097775*01 

>8  8 35  3 | *0 1 

>097775*01 

.5570909*0 | 

•3187505*01 

,5570910*01 

,5990539*01 

,3187505*01 

,5»*0538»0| 

.2180385*01 

,2180385*01 

■677*1 9**0 | 

. 1 107015*01 

,67781 99*01 
.6375001 *01 

, | 107015*01 

•6375000*0| 

.8795855-05 

>795857-05 

floor  SECTION  NUMBER 

y 

II  X*  ,6687063*0? 

J 

BOOT  section  number 

1?  X«  .7105769*02 

j 

,0000000 

L 

.6375000*01 

>000000 

* 

>375000*01 

. 1 107006*0| 

.677*150*01 

| 1 107006*01 

>2781  99*01 

.21*0377*01 

.599059 1 *01 

,71*0377*01 

,5990591*01 

,31*7998*01 

.5520913*01 

>|87997*0| 

>570913*01 

,809776B»0| 

.9*83536*01 

.9097768*01 

,98*3536*01 

.98*3530*01 

.8097775*01 

j 98*3530*0 1 

>097775*01 

,8570*09*01 

.3187505*01 

,5870909*01 

,3187505*01 

.89905 18*0 | 

. 2 1 80385*0 1 

>990538*01 

>1*0385*01 

. 6778 1 98*0 | 

. 1 107015*01 

,6778198*01 

] 1 107015*01 

.6378000*01 

.8798856-05 

>375000*01 

>795855-05 

floor  section  MUmrER 
r 

. onnoooo 

,1107006*01 
,?|80377*0| 
,3I*7997*0| 
.9097768*0j 
. 9*8  3530*0 i 
.S57o’09*O| 
,899o538*0| 
.87781 M 8 * 0 | 
,6375000*0| 


13  X-  ,757997S*0? 
2 

.8375000*01 
, 827R | 89. o I 

.5990581 *01 
.55709|3*0I 
.8883538*01 
. 8097775. oi 
.3187505*01 
.7180385*01 
. I 107015*01 
,8785855-05 


BOOT  SECTION  NUMBER 

r 

.onnoooo 
1 I I 07008*0  I 

>1*0377*01 
>1  *7997*01 
>097768*01 
] 8883530*01 

• 55  20909*0  t 
>99o53fl*OI 
,877R|88*oi 

• 8375000*0  I 


18  X«  .7983181*02 
7 

.8375000*01 

>??8|99*0I 

>99059 | *0  I 
.5570913*01 
.8883538*01 
>097775*01 
] 3 t 87505*0  I 
>180385*01 
. I 107015*01 
>785855-05 


J 


floor  SECTION  NUMBER 
r 

.oonoooo 

. I I 0 7008  *0 | 

. 2 I 8q3  7 7 ♦ 0 | 

, 3 I 87897*0 | 

, 8097788*0j 
,8883530*0] 
.5570’09*0| 
,599o538*0| 

, 8778 | 88*0 | 
,8375000*01 


15  X*  .8381888*02 

7 

.6375000*01 
,6778|89*oi 
.5990581*01 
.552091 3 *0 | 
.8883536*0 j 
.8097775*01 
,3187505*01 

.2180385*01 
. I 107015*01 
.8785855-05 


BOD T SECTION  NUMBER 
y 

>000000 

J I 107006*01 

,2|8o377»OI 

>187997*01 

>097768*01 

,8883530*01 

,5520909*01 

,5990538*01 

>778  I 8**01 

.6175000*01 


16  X*  .8780598*02 
2 

>375000*01 
>2  78  | 99*01 
.599059 |*0| 
>520913*01 
>883536*01 
>097775*01 
,3187505*01 
,2180385*01 
, 1107015*0) 
.8795855-05 


— - ^ 


mis  PAGE  IS  BEST  QUALITY  PRACTICABLE 
HUM  COPY  FURNISHED  TO  DDC  ^ 


Boor  SECTION  NUMBER  17  *•  , 9 I 99300*02 


BOOT  SECTION  NUMBER 


•961 0006*07 


,0000000 

. I I 07006*0 | 

,2l*0377*0j 

,3I87997*0| 

,909776B»0| 

,9**3530*0| 

,557o’09*Oi 

.599053»*0| 

,627B|98*0| 

• 6 3 7 S000  *0 | 


• 6375000*0 1 
,627*199*01 
.S990S9|*oi 
•5520913*01 
. 9**3536*0 | 
.9097775*01 
•3|B7S05*0I 

• 2 I BO  385*0 j 
. I 107015*01 
•B795855-05 


{ooooooo 

, 1)07006*01 
y2l«p377«OI 
,31*7997*01 
{*097768*01 
,9**3530*01 
, 552Q909  *0  I 
,5990538*01 
,627*198*01 
.6375000*01 


,6375000*01 
.6278 | 99*01 
.5990591*01 
,5520913*01 
{ 9883536*0  I 
,9097775*01 
^31*7505*01 
,21803*5*01 
. I I070I5*0| 
.8795855*05 


BOOT  SECTION  NUMBER  |9  X.  ,1003671*03 


BOOT  SECTION  NUMBER  20  *■  .1095592*03 


, 0000000 
. I I 0 7006 *0 | 
,2I*0377*0| 
. 3 I *7997*0 | 
,909776S*0| 
, 9**  35  30*0 | 
,557o’09*0| 
,5'»«n538*0| 
,6?78|98*0| 
, 6375000*0 | 


.6375000*01 
.627B|99*01 
.5990591*01 
,55209|3*01 
, 9**3536*01 
.9097775*01 
. 3 1*7505*01 
, 2 I *03*5*01 

. I 107015*01 
.*795855*05 


jOOODOOO 

. I 107006*01 

, 2) *0377*01 

,31*7997*01 

j*09776**0l 

,9**3530*01 

.5570909*01 

>9905  3 8*01 

{62781*8*01 

16375000*01 


.6375000*01 
, 677«|99*0| 
J 5 9905  9 | *01 
.5570913*01 
,9*83536*01 
. 9097775*0  I 
.3187505*01 
.71*0385*01 
1107015*01 
>795855-05 


BOOT  SECTION  NUMBER  2|  X*  .1087913*03 


boot  section  number 


, I I 79283*03 


, 0000000 
. I I07006*0| 
,21*0377*0| 
.31*7997*01 
,909776**0| 

, 9*8  35  30*0  I 
.5570909*01 
•599o53»«0| 
.677*198*01 
•6375000 *0i 


.6375000*01 
•627B|99*oi 
.5990591*0) 
.55209 | 3*0 1 
.9**3536*01 
.9097775*01 
.3187505*01 
.2180385*01 
. I 107015*01 
.*795855-05 


{ooooooo 

, I (07006*01 
,21*0377*01 
,31*7997*01 
,9097768*01 
{ 9*8  35  3o  *0  I 
,5520909*01 
{5990538*01 
.6278198*01 
« 6 375000*0  I 


{6375000*01 
,6278 |99*0I 
,5990591*01 
>52091  3*01 
>**3536*01 
|9097775.0l 
.3187505*01 
{2180385*01 
. I 107015*01 
>7*S855-05 


Bony  SECTION  NUmbER  23  X.  .11711 59*03 


BOOT  SECTION  NUMBER  29  X*  .1209277*03 


.0000000 
. I I 07006*0 | 
.21*0377*01 
,3I»7997*0| 
, 909  7768*0 | 
. 98*3530*01 
•55?o’09*0| 
, 59*0538  *0 | 
,627*|98*0| 
,637SOOO»0| 


,6375000*01 

, 62781*9. oi 

.5990591*01 
.5520913*01 
.9883536*01 
. 9097775. oi 
.3|»7505*0| 
.21803*5*01 
• I I070|5*0I 
.8795855-05 


jonooooo 

{ I 107106*01 
,21*0573*01 
,3|*7785*0| 
.909*138*01 
>8*397  j*Ol 
{5521*07*01 
,599107**0) 
,6278719*01 
{6375575*01 


« 6375575*0  I 
{6278716*01 
,599|0«l*0| 
{552191 |*0| 
,9883976*01 
{*098199*01 
,31*7792*01 
,2180581*01 
,1107(15*01 
>796699-05 


r • 


IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
raw*  COPY  JTJRtilijHED  TO  DLC 


Boor  SFCTION  NUmbFR  25  X*  ,1212434*03 


BOOT  SECTION  NUMBER  24  X. 


. | 27 | 228*03 


y 


I 


Y 


2 


.0000000 

. I 0*1 I 5 l *0  I 

• 2 I 55057  *0 | 

• 3 I 50*8  3*0 i 
,10501 *3*0 | 
, l8?487o*0| 
.5154777*01 
,5720,73*0| 
.4705273*01 
, 4 300,70*0 | 


•4300770*01 
•42052l5*oj 
.5720774*01 
.5154802*01 
.1824824*01 
.1050187*01 
.3150170*01 
.2155045*01 
. i071|40*0l 
. 8411271-05 


iono 0000 

J 1047221*01 

,2105741*01 

,3078708*01 

.3757712*01 

i 17 i 4854*0 1 

|53321*i *01 

j57»4082*0l 

,4043875*01 

.4157122*01 


,4157122*01 
j 404  38  7 7 *0  I 
,5784087*01 
^5332185*01 
,1714841*01 
^3157718*01 
^3078714*01 
j 7 I 05748*0  I 
, 1047237*01 
.8117340-05 


Bony  SFCTION  NUmbEP  27  X.  .1275054*03 


BOOT  SFCTION  NUMBER  28  X. 


.131111 7*03 


r 

, 0000000 

. 1 01 1 217*0  I 
. 205084  1 *0  I 
« 777R 1 58*0 1 
,3*51 358*0 | 
,151J115*0| 
,5172743*01 
,543i4i4*0| 
.5105227*01 
,5l7»37l*0| 


Boor  SECTION  NUMBER 


,0000000 

,177o737»Oo 

, | »21157*0 1 
,2*1 337o*0| 
.3*14777*01 
,l3lo33l*Oj 
.1*72701*01 
,52*7107*01 
,55i|24|*0| 
.5*74715*01 


2 

.5774321*0' 
.5705223*01 
.5431477*01 
.5177747*01 
.1573150*01 
.3851341*01 
.7778145*01 
.7050848*01 
. 1011257. oi 
.8274314-05 


77  X.  .1328111*03 
2 

.5424715*01 
.551)242*01 
.528711 | *01 
.1*72705*01 
.1310331*0| 
•34)4805*01 
.7813377*01 
. j727l44*0| 
.7770807*00 
.7717325-05 


BOOT  SECTION  NUMBER 

y 

,0000000 
.737434  I *Oo 
, |850722*0| 
,2705575*01 
,317827o»0| 

. 1 1 15  I * I *0 1 

,14B4|11*0| 

1 .50*1817*01 

.5378115*01 
,51|  I 1 5j  *0 | 

> 


t 


31  X.  .1312710*03 
2 

.5111 153*01 
,5328714*0i 
.5081821*01 
.1484177*01 

.1H5|84*0I 
.3178225*01 
.2705581*01 
. 185072’*0| 
.1374131*00 
.7173555-05 


y 

0000006 

1008054*01 

1785183*01 

7*02587*01 

373|1*8»OI 

1117015*01 

5027121*01 

5155072*01 

5714171*01 

5*05141*01 


2 

,5805141*01 
J57 14174*01 
,5855075*01 
i5027l25*01 
jl*170 20*01 
,373|111*01 
,2102587*01 
, 1 78511 1*01 
j 1008041*01 
.7741105-05 


BOOY  SECTION  NUMBER  30  *• 


. | 337HS*03 


y 


2 


0000000 
7528175*00 
1874488*01 
2713531*01 
3577070*01 
1703337*01 
175|717*01 
5 I 54  I 4p*0 I 
510371 |*0| 
,5187073*01 


i5187073*01 
j5*037 I 2*0  I 
,5154143*01 
j 175 1 114*0 1 
,1203311*01 
j 3527024*0  I 
,2713511*01 
, 1874475*01 
,1528248*00 
,7527710-05 


BODY  SECTION  NUMBER  32  *« 


*1317174*03 


y 


2 


3401317*00 
1200*53*01 
20107*4*01 
774147 j *01 
31582*5*01 
*057025*01 
1554740*01 
1138811*01 
5 | 151*7*01 
5370247*01 


,5320241*01 
^ 5 1 1511 1 *0  I 
,1138817*01 
1 *554747*01 
,1057027*01 
,3158250*01 
.2741*77*01 
j 20 1 0717*0 1 
,1700840*01 
« 3401 3*7 *00 


j 

1 


42 


THIS  PAGE 

FROM  COP* 


,SBfSTQU«lI”mCIiaB“ 

FURNISHES)  TO  M>C  


boot  section  number 

Y 

,*T727*0*Oo 

, | 27*999*0  | 
.20*7*78*01 
,27*| 32**0| 
.3*1**I2»0| 

, 3***7A3*0| 
,**4*|B**0| 
,*S3553|*0| 
,50*1**0*01 
,5227*3**0 | 


Boor  SECTION  NUMBER 
T 

, 47*7828*00 
, I 3 | * 3 I 0*0 [ 

, | *8  3*3A*0 | 
,2*0722**0i 
,3l75237*0| 

, 3*75805 *0| 

, *0*82*8*0 | 

, *833**2*0 | 
,**75013*0i 
, *8 | 4*85*0 | 


BOOT  SECTION  NUMBER 
T 

,5»5*2***Oo 
. | 1 8**5**0 | 

, | 788*7**0 1 
,23**378*01 
,28*0*78*01 
. 33l0**l*0| 

, 3** I 2 I 8*0] 

, 3’*32*B*0  j 
, * 2 1 0*  28*0 1 
.*338751*01 


BOOT  SECTION  NUMBER 
Y 

, *078007*00 
.**'51 30*0o 
,15023*7*01 
.20070*8*01 
,2***2I5«0| 
,2***3*3»Oi 
,3207  358*0 i 
.3*7752**01 
,3*8b8*|*0| 
,37*707**0i 


33  X.  .1357753*03 
2 

.5227*38*0  I 
.80*1*87*01 
.*835538*01 
.****188.01 
.3*8*7*8*01 
.3*l**l«»01 
.77*1132*01 
.70*2*38*01 
. 1275005*01 
,*7727*0*00 


35  X.  .1371*83*03 
Z 

, *81*985*01 
.**750l'**OI 
.8*337*5*01 
.*0*8301»0| 

. 3*7580**0 1 
.31 752*2*0i 
.7*07231*01 

.1*838*2*01 
. I 3|8li 8*0  I 
.42*782**00 


37  *«  .13*0**5*03 

Z 

.*338751*01 
.*210430*01 
.3**3251*01 
.34*1222*01 
.3310**5*01 
.2840*82*01 
.73**382*01 
. j 788*80*0) 

• 1 16*4A**0| 
.545*2*8*00 


3*  X.  .1*1 0705*03 
Z 

. 37*207**0 I 
. 345*8*7*01 
.3*77574*01 
•3207340*01 
.784*3*7*01 
.2*4421**01 
.7007073*01 
.1807357*01 
.9*38180*00 
.8028007*00 


BOOT  SECTION  NUMBER 
Y 

>90**58*00 

. I 3 79  | 8 * * O | 
.7018100*01 
, 7 702  I 2 | *0  I 
, 370*473*01 
J 38182*1 *01 
. *2"5*15«01 
>438 1*4*0) 
,*888713*01 
.5031*00*01 


BOOT  SECTION  NUMBER 
Y 

j* 1 17277*00 
, I 748* I5*0| 
, 18*88*7*01 
j 2*8* I 8**0  I 
,3027220*01 
,3801*52*01 
,390251 1*0) 
,*221 37**OI 
,**518*5*01 
>588185*01 


BOOT  SECTION  NUMBER 
Y 

;***I*13*00 
j 1087102*01 
, 1*54*02*01 
,2189451*01 
.7*752*3*01 
>102408*01 
,3*42**4*01 
,37*4978*01 
,3**9851*01 
.8(14*81  1*01 


BODY  SECTION  NUMBER 
Y 

, ?*3*A2**00 
>27*201  *00 
j | 3.1 1 31  i *01 
jl8080*i*0l 
,77*10*8*01 
,7470009*01 
,7935433*01 
, 11*0 3 *0*01 
i 1 1*B2  2 1 *0 l 
>*1521  9*0  i 


3*  X • .13418*3*03 
Z 

,5031400*01 
,*888715*01 
, *438 | *8*0 | 
j 1285* I * » 0 i 
, 38382*4*01 
| 3304426*01 
,2702| 74*01 
j 2038 1 04*0  I 
.132*1*3*01 
. 5*0**58*00 


34  X.  .1381182*03 
Z 

j*688 (85*01 
j *85 15*7*01 
,8721378*01 
,3*02515*01 
,3501454*01 
,3027225*01 
^ 2*8* | 8**0 | 
,18*8853*01 
1 248421  *01 
> 1 1 7277*00 


38  *■  .1*00782*03 
Z 

,804481 1*01 
>9**852*01 
, 37*4*24*01 
, 3*A2***»01 
,3107417*01 
,2*752*7*01 
,718**54*01 
>454*07*0) 
>087107*01 
.8**1  * I 3*00 


*0  X.  .1*2077**03 
Z 

>*3521**01 
>3*8727*01 
>1803*7*01 
>935435*01 
>42001  2 • 0 I 
>2*I0’7*01 
, 1 8080*5*01 
>331315*01 
.872*7*7*00 
>93442**00 


IBIS  PAGE  IS  BEST:  QUALITY  PRACHCAilLl 
oofY  PURUilSHHD  TO  DOC  — " 


- 


BOOT  SECTION  NUmbFR 

*1  *•  *1 *308*3*03 

BOOT  SECTION  NUMBER  *2 

X*  .1**104**03 

Y 

. 1 4 1 7 1 2 1 *00 

2 

•3101084*01 

Y 

17*3*229-07 

L 

• 27  409*0*0 1 

, 447?8"**Oo 

• 30338  24*0 i 

1*7*3537  *00 

.7718882*01 

.1  1*58I0*0| 

•28841*7*01 

j9»*|*2**00 

| 28**01 7*01 

,|5*8*59‘0| 

• 744  1 94*  *0  | 

I 1 3*02*9*01 

,73*0488*01 

, 200*725*0) 

.7347220*0) 

1 177**07*01 

,21 1 *489*0 | 

,2347218*0| 

.7009728*01 

,71 1*487*01 

, 177**10*01 

, 244 | *42*0 1 

. j598943*0l 

j 73*048 J*OI 

,13807*7*01 

.2«*4I*S»0| 

.1 1 *88 1 * ♦ 0 i 

,25**01 )*0I 

,***1*89*00 

,J0J|B2S*0| 

•44229 30 *0Q 

,771 888 | *0 | 

,97*3873*00 

.3101 0*4*0 | 

•1412121 *00 

•7740**0*01 

, 7*3*22*-07 

BOOT  SECTION  NUmbFR 

v 

*3  X«  • 1 9*70* ^403 

body  section  number  ** 

tm  .1*53029*03 

* 

T 

,0000000 

2 

.2558341*01 

Y 

loOOOOOO 

l 

j 2354  3 1 2*0 1 

,***75*8*00 

.251**94*01 

• *0* 1 4*0*00 

17320515*01 

• B 7 So  1 1 1 *00 

.7*0*074*01 

j 805*054  *00 

,221*209*01 

, 1 77*| 8| *0| 

.2715408*01 

, 1 1 78 |55*0I 

,70*0427*01 

, 14***83*01 

. |*5*87I*0I 

j 151*407*01 

1 1 8050* | *0 1 

, 1 *5*8 1 9*0  j 

. 1 4*9*84*01 

, 1 *0503**0 1 

,151*410*01 

•27I5407*0| 

. | 27  9 | 8*  *0  1 

1 70*0475*01 

,1 178158.01 

,?*0*07m*0| 

.87501*9*00 

^771*208*01 

1805*0*4*00 

.251***4*01 

.9**7582*00 

. 7370519*01 

.*0*1772*00 

•2558343*0i 

■ 350*8 1 4*05 

•7354312*01 

13232427-05 

BODY  SECTION  NUMBER 

*5  X*  .1*5*017*03 

body  SECTION  NUMBER  *4 

X.  .|*4****» 03 

Y 

7 

Y 

7 

,0000000 

.7I5*2**«0I 

lonooooo 

1 1*52217.01 

.37*08*0*00 

.2121544*01 

1 338**87*00 

1 1*22559.01 

,73481 1 **Oo 

.207*375*01 
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The  FORTRAN  listings  are  preceded  by  a list  of  the  nonsystem 
subroutines  required  by  each  of  the  three  main  programs 

Program  1: 

Main  Program  AXISYM 

Subroutines  AXI3D 

BDYAX 

LAGINT 

PLOTA 

Program  2: 


Main  Program 

APNDG1 

Subroutines 

APPINT 

INTSEC 

APPNDG 

LAGINT 

BODY 

LINE 

DELBDY 

PLANE 

DOUGC 

FINCRD 

F INEND 

PLOTC 

Program  3: 

Main  Program 

OFFBDY 

Subroutines  Nonrequired 


IBIS  PA&g  IS  BEST  QUALITY  PRACXICABU 
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C 

C PROGRAM  I - AXlSYM 

C 

C main  PROGRAM  - AXISYMMETRlC  BODY  WITHOUT  appendages 

c 

DIMENSION  LABEUHOI 

DIMENSION  XB0DY(10O,1O>,yB0DY(1OO,IO),ZBODY(|O0,IO) 

DIMENSION  X ( 200 ) ,Y(200> 

N I »5 
N0«6 

nP-1 

I z • X 0 

READ ( N I ,50  I ) LABEL 

rEAD(N1,5O2)JREAD,IPRINT,IPL0T, I PUNCH, NS YM, IDOUG,NPtS,NBODY 

wR1T£(N0,530)NSYM 

WrITe(N0,S03)LABEL 

wRITe(N0,50M) 

wRITeINO.SOS) IREAD,IPRINT.IPUNCH,NPTS,NP0DY,IPL0T,ID0UG 
lP(NsYM.GT,3,0R.N5YM.LT,i)WRITE(N0,53ll 
IP( Np TS. GT, 200 > WRITE! NO, 5O6) 
lP(NpTS.LT.20)WRITE(N0,533) 

IP ( NbODY ,GT . 100  >WRI TE  » NO  ,S07  ) 

IF(NsYM,gT.3.0R.NSYM.lT,i)ST0P 

IP ( NpTS.GT, 200. OR, NBODY, gT, 100 (STOP 

lP(NpTS,LT.20)ST0P 

GO  To  ( 1 , 2 ) , I READ 

1 DO  3 1*1 ,NPTS,3 

3 READ i N! ,508 >X I I ) , Y ( I ) ,X I l*\  ) , Y(  I ♦ I ) ,X ( 1*2  1 , Y(  1*2  ) 

GO  To  S 

2 READ(NI,508)(X(I),I*I,NPTS) 
rE  A D ( N I ,508)  ( Y I I ) , I a 1 ,NPTS) 

5 CONTINUE 

WR I Tf ( NO , SOP  ) 
wRITe<N0,510) 

DO  6 1*1  , NP  T S 

6 wRITe<N0,S)  l ) I ,Xt  I ) ,Y(  I ) 

CALL  BOYAXIX, Y,NPTS,XB0DYiYB0Dy,ZB0DY, NBODY, IZ.NSYM) 

I P ( IpRINT.EG.01G0  TO  1 1 

wpiteino.sis) 

DO  1 6 I»1  , NBODY,  2 

I I » l ♦ 1 

I P ( I I • GT .NBODY ) GO  TO  lS 

wRITe<NO, 516)1, X BoO  Y ( I , X ) . II, XBODY ( I ♦ 1 , 1 > 
wRIT£(N0,S17) 

00  I 3 J»  1 , I Z 

1 3 WpITE<N0,5l8)2B0DY(I,J),YB0DY(I,J),ZB0DYtI*l,J),YB0DY(Ul,vn 

GO  To  16 

JS  WR I Te ( NO ,51 9 ) I ,XBODY(  I , I ) 

WrITeIN0,520) 

DO  17  J* 1 , I Z 

17  WRITEIN0,521)YB0DY(I,JI,ZB0DY(I,J) 

16  CONTINUE 

11  CONTINUE 

1 F ( I DOUG. EG .0 ) GO  TO  90 
wRITe<NO,SOO)NSYM 

wRITe(N0,S03)LABEL 
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C * L L AX  I 30( XBODY , YBoDY ,ZflOOY ,NBOOY , 1Z , !PUNCH,NO,NP ) 

*40  I r < IpLOT.EQ.OIGO  TP  50 

W*  I T|r  ( NO  , 502  ) 

CALL  PLOT  A ( XBODY , VBOOY  , zboo  snbody  , 1 Z ) 
so  continue: 

STOP 

boo  forma t ( i h i .sx , * coordinates  for  input  to  the  Douglas  three  dimen 

ISIONaL  POTENTIAL  flow  PROGRAM'  ,/,  5X  , ' HAVING  * , I •*  , 2X  , ' PLANES  OF 
1SYMMETRY'  ,///  ) 

50  l fORMAT  ( <i  0 A 2 ) 

502  FORMaT(8I3) 

503  FORM  a T ( 2X  , *40A2  , / // • 

5 0 *4  FORMAT  ( 1 MX  , ' CONTROL  PARAMETERS'  ,/// I 

505  FORMAT ( zt  ,'  IREAD  = ' , I H , 2 X , ' I PR  I N T « » , J H , 2X  , * I PUNCH  * 

13xi'NPTS  = ' , IM  ,3X  ,'NBOOY  ■ ' , I *»  , 3 X , ' I PLOT  » ' , I H i 3 X , * I DOUG  » 

I I **•///  1 

506  fORMAT I 2X  , ' NUMBFP  OF  INPUT  POINTS  EXCEEDS  200  - PROGRAM  TERMINATE 

lEO*  ,//  ) 

507  FORMAT ( 2X  NUMBER  OF  OUTPUT  POINTS  EXCEEDS  100  - PROGRAM  TERMJnAT 
lEO'  ,//  ) 

SOS  FORMAT ( 6F 1 0.6 ) 

509  FORMAT ( 1H  1 ,9X  , 'INPUT  BODY  COORDINATES',//) 

5 1 C FORMAT ( i MX  X'  , 1 6X  /'  ,// > 

511  fOPMa T ( 2X , l 3 , 2 ( 2X  ,F  l R - 7 ) | 

515  fORMa T ( i h 1 , 5X CALCULATED  THREE  DIMENSIONAL  COORDINATES',//) 

516  FORMAT  (/// ,2X  ,'  BODY  SECTION  NUMBER  ' , I R , 2 X , ' X « • , £ 1 *»  . 7 . *♦  X , ' BOD  y SECT 

1 I on  NUMBER'  , 19, 2X  ,'XS'  » E 1 B , 7 ,/) 

517  FORMlTl  iuk.'y'  ,17X('Z'  , 2fl>!  , ' Y ' , 1 7X  , ’ Z'  , / I 

513  FORMaT(2(MX,L1‘4.7),12X,2(BxiEih,7)) 

519  fORMaT  (///,2X, ‘BODY  SECTION  NUMBER  ' , I *4 , 2 X , ' X » » , E | H . 7 i / I 
5 2 0 FORMAT ( 1 QX  ,' Z'  , 1 7X  ,' Y'  ,/ ) 

5 2 1 FORmAT(2IRX,E1<4.7)) 

530  FORMa T ( 1 H i , 4X , ' C ALCULAT I oN  OF  THE  THReF  DIMENSIONAL  COORDINATES',/ 
1,5X,'F0R  AN  AXISYMMETRIC  »ODY  WITHOUT  A P P E ND A GE S ' , / , 5 X , ' H A V I N G * , 

2 I R , 2 X , * PL  A N£  5 OF  SYMMETRY'  ,///) 

5 3 1 FORMa T ( 5X NUMBER  OF  PLANES  OF  SYMMETRY  INCORRECT',//) 

532  FORMAT t 1 H 1 ,5X ,' PLOT  INFORMATION*  ,//) 

533  FORMAT ( 2X  , * NUMBEP  OF  INPUT  POINTS  IS  LESS  THAN  20  - PROGRAM  TE*H 
1 jNATeD'  ,//) 

p Nf) 
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PROGRAM  2 . APNDG1 

MAIN  PROGRAM  - AXlSYMMETRlC  BODY  WITH  APPENDAGES 
TWO  0R  THREE  PLANES  OF  SYMMETRY  ALLOWED 
dimension  label ( MO  t 

DIMENSION  COEF!M) 

DIMENSION  XB | 25  ) , HC ( 25 ) 

OIMEnSION  X ( 200  I , Y ( 200  I 

DIMENSION  XAPP! I M , 1M  > , YApP<  I M , i M I ,ZAPp!  | M , JM  ) 

DIMENSION  XINC(  1 M ) ',  Y 1 N C ( | M ) , Z l N C ( 1 M ) 

DIMENSION  XBODY(  100,10  1 ,Y0ODY(  100,10)  ,ZROOY( 100,101 

N I »5 

no»6 

nP»1 

I I 0 

READ) NI  ,201  (LABEL 

rEAD(NI,202)IREAD,IPR1NT , I PLOT , IPUNCH.NSYM, I DOUG , NBoO Y , NF I N , JFWD, 
1 I AFT  * I F I N 
wRITe(N0,221  (NSYM 
W«ITeIN0,217 (LABEL 
WP I Te ( NO , 2 1 8 ) 

wPITe(N0,2|6(IREAD,IPRINT.IPL0T. IPUNCh,N80DY,NFIN, IFWD,  iaft,nsym 
1 »I00UG,IFIN 

lF|NSYM,GT.3.0R.NSYM.LT,2lWRITEIN0,2l9) 
IF(Nb0DY.GT.200>WRITEIN0,300( 
lF< NbODY. L T.20> WRITE« NO ,352 ( 
lF(NplN,GT.25)WRITE(N0,30l » 

I F ( IfIN.GT, |N( WRITE! NO »350) 

IF(NeIN,lT,5»WRITE<N0,353( 

I F ( Ip  IN.LT ,5  > WR I TE  < NO  , 3S  | > 

lF(NsYM,GT,3,CR»N5YM.LT.2»ST0P 

IF ( IdOUG.EQ.OIGO  TO  MO 

G°  To<36,36,37> ,NSYM 

IF ( lEWD.GT.M0)WRITE(N0,305» 

l F ( IaFT.GT.IO)  WRITE<N0,307» 

JF ( IpWD.GT.MO.OR. I AFT.GT, 10 (STOP 
GO  To  38 

lF<Ipw0.GT.30>WRITE(N0,302> 
lF(  IaFT,GT,20» WRITE! NO,  303» 
lF ( I pwD ,GT , 30  # OR  « I AFT, GT, 20 (STOP 
IF  < IFWD.GT ,60 (WR ITE ! NO, 30®  > 

IF | lAFT,GT,2b)WRITE!N0,309» 

IF ( lFWD.GT.60.0R, IAFT.GT.2BIST0P 
CONTjNUE 

IF(Nb0DY,GT.200.0R.NFIN,gT,25(ST0P 

JF ( Nfl ODY.LT. 20. OR, NFIN.LT. 5 ( STO9 

1F(1FIN.GT.IM.0R.1F1N.LT.5)ST0P 

LPTS.IFIN 

L5ECT" I F I N»3 

GO  TOI205.206)  , 1READ 

DO  1 I » I , NBOD Y , 3 

RE  AO ( N I , 2 > X ( I ) , Y I I ) , X I I ♦ | ) , Y I I ♦ 1 I , X ( I«2) , Y ! 1*2  1 
GO  To  207 
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206  READ(NI  ,2 M X(  I I , l«l  ,NB00Y  I 
READ | N I ,2  ) < Y 1 J ) ,1*1  , NBOD Y > 

207  cONTjNUE 
wRITe(N0,3) 

WR I Te ( NO , 22  ) 

DO  1 9 1*| , NBJDY 
19  WR 1 Te ( NO  , 4 ) I , X I I ) , Y ( I ) 

wR I Te ( NO  , 1 Ml 

9EAD(NIf<))CORDT,COPDR,XSTART,XOESFT,YTlP,TCREE,TCTlP,TCROOT 

look.o 

DO  ? 2 0 1*1, N BODY 

IF(LoOK,NE.O)GO  TO  22 0 
220  jF ( X | I > ,GE . XST ART  )LOOK* I -1 

CALL  LAGINT(X,NB00Y,Y,R,xSTART,YR00T,L00K,lS,C0EF) 
wRJTe:N0,1 1 IXSTART.XOESETiYTIP.CORDT.COPOR 
wRITe(NO,12)TCREF,TCTiP,tCROOT 
GO  To ( 2 1 1 I 2 1 2 If  I RE  ad 

211  DO  2 i 3 1*1  ,NF1N,3 

213  REaD(NI,2>XB(I>,HC<I>,XB(U1),HC(I*I>,XB(I*21,HC<I*2» 

GO  T o ? I H ' 

212  rFaD(NI  ,2)  (Xb(  I ) , 1*1  , N F I N ) 

rFAD|NI  ,2 ) (HC ( I ) i I« 1 ,nfin» 

21H  CONtINUF 

WR I Te ( NO , 24  ) 

DO  2 J I * 1 , NF  1 N 

21  wRITeINO,  13)  I , X B ( I ) ,HCU  ) 

CALL  APPNOGl X , Y jNBODY ,XB ,HC  ,NF in ,XST Art ,XENO  .XOFSET  , Yt IP  , YROOT  , 
1C0RDTiCORDR,L5ECT ,LPTS,XINC,YINC,2INC,XAPP,YAPP,ZAPP,MSECT, 

2 T 0 RE  r , TCT I P , TCRO0T  1 

43  CALL  BODY(X,Y,NBODY,XAPP,YAPP,ZAPP,K,lZ,MSECT,LPTS,XBoDY,YBODY, 
1ZP0Dy,XSTART,IFWD,IAFT,NSYM) 

44  c ON  T I NUE 

I F ( IpRINT.EQ.OlGO  TO  1 10 
WR 1 Te ( NO  , 82  ) 

WRITE (NO,  88) 

DO  83  L* 1 ,lPTS 

03  WR  1 Te  ( NO  , 3E)  )L  . X I NC  ( L ) , Z INC  { L I , y INC!  L 1 
wR I Te ( NO  , 1 5 ) 

DO  4 0 0 ! = 1 ,MSECT  , ? 

I I * I ♦ 1 

IF ( I J ,GT . H S F C T ) GO  TO  401 

WRITeINO,  18)1,  YAPP(I,l),!l  .YAPPIUI.i) 

WR I TE ( NO , 25 ) 

00  2 9 J* 1 ,LP  rs 

29  WRITeINO,  16)  XAPP1  I ,JI  « Z APP ( I ,J)  ,XAPPj  I ♦ | , J ) ,ZAPP1  I ♦ I , J ) 

GO  To  400 

4 0 I WR I Te ( NO , 402  i I , YAPP ( I , 1 ) 
wR I Te ( NO , 403 ) 

DO  4 0 4 J * 1 , L P T 5 

404  wRITe<N0,405)XAPP(I,J).ZaPP(!,J| 

400  CONT  j NUE 

WRITEINO, 5) 

DO  loO  I = l , K , 2 

1 I * I ♦ 1 

I F ( I 1 . GT.K 1 GO  TO  | 0 i 

WRITE (NO ,20)  I ,XBOOY(  I , I > , I I , X B OD Y ( I ♦ j , 1 ) 

WRITE (NO, 23) 
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do  i o o j»i  , iz 

WRITER  NO,  8)  ZBOOY ( I , j) .YBODYI  l , J) ,ZBODY( I ♦ | , J ) ,YBOOY(  1*1  , J • 
60  To  102 

WR I TE ( NO , 30  > I ,XBODY(  ! , I) 

WRITE(N0,31  ) 

DO  3 2 J ■ 1 . I Z 

wRITe(NO,33)ZBODY( I , J > .TODDY!  I , J ) 

CONTjNUE 

lF(  IpL0T,EQ.0)60  TO  130 
WRITEINO.M  ) 

CALL  PLOTCIXBODY.YBODY.ZBODY.K, 1 Z , X A PP s Y A PP , Z APP , MSE C T , LPT S ) 
CONTjNUE 

I F ( IOOUQ.EQ,0)60  TO  140 


WRITE  I NO, 2 15) 
wRITe(N0,217)LABEL 

CALL  D0U6C  ( XBOD Y , Y BOD Y , Z BOD Y , K , I Z , X APP , Y A PP , Z APP , MSEC T , L P T S , 

1  I PUNCH  , NO  , NP  ) 

1 <40  CONTjNUE 

2 P0RMaT(6F10.6) 

3 FORMAT ( 1 H 1 ,9X  ,'  INPUT  BODY  COORDINATES',//) 

4 FORMAT  I 2X , I 3 , 2 < 2X  ,E  1 4 .7  ) ) 

5 FORMAT< | H 1 ,5X, 'CALCULATED  BODY  COORDINATES',//) 

6 FORMAT  I 8F I 0.6 ) 

0 F0RMaT(2(MX,EM.7),I2X,2(*)X,E14.7)) 

9 FORMAT!//) 

11  FORM4T(2X,'XSTART»',E1,).7.2X,'XOFSET»',F1‘).7,2X,'ZTIP.',E1*).7»//i3 
IXi'CoRDT-'  .EIR^.SX.'CORDR**  ,E  1*4.7,/) 

12  F0RM4T(3X,'TCREF«'  ,EI4.7,3X,'TCTIP»',E14.7,3X,'TCR00T»',E14.7,//) 

13  FORMAT  ( 2X  , I 3 , 2 < ^ X ,E  1 *4,7  ) ) 

14  FORMAT!  1H)  , I2X, ' INPUT  APPENDAGE  DATA',///) 

15  FORMAT( 1H1 ,2X, 'CALCULATED  APPENDA6E  COORDINATES',//) 

16  f0RMaT!2!2X,E1‘),7),2OX,2|2X,E14.7)) 

18  FORMAT!///, 2X/APPENDA6E  SECTION  NUMBE R ‘ , I 4 , 2 X , ' Z - ' , E 1 4 . 7 , 4 X , 

1  'ApPENDA6E  SECTION  NUMBER ' , I 4 , 2X  , ' Z« ' , E 1 4 . 7 , / ) 

20  FORMAT!///, 2X, 'BODY  SECTION  NUMBER ' , I 4 , 2 X , ' X . ' , E 1 4 , 7 , 4 X , * BOO Y SECT 
1 I ON  NUMBER'  ,I4,2X,‘X«'  iEl4,7,/> 

22  FORMAT ( 1 4X  ,' X'  , 1 6X  ,' Z'  »// I 

23  FORMAT! 10X,'Y',17X,'Z',2RX,‘Y',17X,'Z',/) 

24  FORMAT ( J5X  ,' XB'  , 1 6X  HC'  »// ) 

25  FORMaT(8X,'X  ' » 1 4 X , ' Y ',3MX,'X  ',14X,'Y  ',/) 

30  FORMAT!///, 2X, 'BODY  SECTION  NUMBER ' , I 4 , 2 X , ' X ■ ' , E 1 4 , 7 , / ) 

3)  FORMAT!  10X,'Y'  ,17X,'Z'  ,/) 

33  F0RMaT12(4X,E14.7)  ) 

35  FORMAT ( 5X  , I 3 , 3 ( 3X  ,E  1 4 . 7 | , / ) 

41  FORMAT ! 1 H 1 ,5X  ,' PLOT  INFORMATION',//) 

02  FORMAT! 1HI ,9X,' ITERSECTION  POINT  SUMMARY  - APPEND A6E  AND  BODY',//) 
88  F0RMaT(18X,'x‘,16X,'Y‘,UX,'Z',//) 

201  F ORmAT ( 4 0 A 2 ) 

202  FORMAT ( HI3) 

215  FORMAT! 1H1 ,5X, 'COORDINATES  FOR  INPUT  TO  THE  D0U6L AS  THREE  DlMENSlO 

1 1N*L  potential  FLOW  PR06RAM'  , /// | 

216  F0RmAT(2X,'  IREAD  • ' , I 2 , 4 X , ' I PR  I NT  « ' , I 2 , 2X  , ' I PLOT  ■ ',12, 2*, 

, 1 'IPuNCh  « ' ,I2,///,2X,'nB0DY  . ' , I 3 , 5 X , ' NF  I N ■ ' , I 3 , 2 X , ' I FWO  ■ ', 

, 2 I 3 , 3 X , ' I AFT  . ' ,I3,///,3X,'NSYM  ■ ' , I 3 , 4 X , ' I D0U6  ■ '.I3.2X, 

3  ' I F j N * ',13,///) 
j 217  FORMAT! 2X ,40A2,// ) 

218  FORMaT (//,  15X ,' CONTROL  PARAMETERS',//) 
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219  FORMAT ( 5X ,* NUMBER  OF  PLANES  OF  SYMMETRY  INCORRECT*,//) 

22\  FORMAT!  )H1  ,5X, 'CALCULATION  OF  THE  T HR  £ E DIMENSIONAL  COORDINATES*,/ 

i,sx,*  for  an  axisymmetRic  body  with  appendages*  ,/,6x, 'having* , 

2 I 4 ,2X ,* PLANES  OF  SYMMETRY*  ,///) 

300  FORMAT!//, 2X,*NB0DY  EXCee^ES  200  POINTS  - PROGRAM  TERMINATED*) 

301  FORMaT (// ,2X NF I N EXCEEDES  25  POINTS  - PROGRAM  TERMINATED*) 

302  FORMaT (// ,2X ,* IFWO  EXCEEDES  30  POINTS  - PROGRAM  TERMINATED*) 

303  FORMaT (//, 2X  ,*  I AFT  EXCEEDES  20  POINTS  - PROGRAM  TERMINATED*) 

305  FORMaT (// ,2X ,* IFWD  EXCEEDES  HO  POINTS  - PROGRAM  TERMINATED*) 

307  FORMaT!//, 2X,*IAFT  EXCEEDES  10  POINTS  - PROGRAM  TERMINATED*) 

308  FORMaT (// ,2X  ,»  IFWD  EXCEEDES  60  POINTS  - PROGRAM  TERMINATED*) 

309  FORMaT (// ,2X I AFT  EXCEEDES  25  POINTS  - PROGRAM  TERMINATED*) 

350  FORMaT ( // ,2X  , * I F IN  EXCEEDES  IH  POINTS  - PROGRAM  TERMINATED*) 

351  FORMaT!//, 2X, * IFIN  IS  LESS  THAN  5 - PROGRAM  TERMINATED*) 

352  FORMaT (// ,2X ,' NBODY  IS  LESS  THAN  20  - PROGRAM  TERMINATED*) 

353  FORMaT! //,2X,»NFIN  IS  LESS  THAN  5 - PROGRAM  TERMINATED') 

402  FORMaT (///, 2X ,' APPENDAGE  SECTION  NUMBER  * , I 4 , 2X  , * Z« * , E 1 4 , 7 , / ) 

403  fORMaT!8X,*X  * , 1 4 X , * Y ',/) 

<405  FORMaT  ( 2 ( 2X  ,E  1 4,7  ) > 

310  STOP 

E N D 


F 80  fc  K A h - CFFbDY 

CALCOLATlOf.  LF  IhE  CFF  BODY  POINTS  FOR  THE  DCUbLAS  THREE 
ulhENSlONAL  potential  FLOW  program 

OlhENS ION  LAlELUC) 

C1HENSI0N  TTIi:01),Z2(1C01) ,XX( 1001 >t VAC 10C1), THETA (1001) 

IPOS=0 

5 

N0  = t 

KP=1 

A =0 

I sl  a=c 

I LAST=i 

H E A p ( ft  1 « 2 G 3 ) LA  i.  E L 
.Nl TE (NO  ,203) 

.MITE  (NO  ,iH  * )L  Abe  L 

k EA  0 CM,  !>0C  > APT  Si,  N ThETA,  X,Z$TART,DFLYZ,TS  TART,  BELT 
If  (NPTSZ  .EC  .LKO  TO  901 
iPOS-I F 0 S ♦ 1 
fcRiTt(NC,2L0)lFCi 

■ Ri Tt (NO.207  )NFTOZ,NTHE  TA ,X,ZSTaRT , DELVZ  tTSTART,DELT 
YR(1)=2START 

THtTA<T)=TSTART*i.1Al5V/1i<0. 

THt TA (1 )=0 
0 0 2 1=2  »NP  T SZ 
YR(  I ) = Y R U-1)4wclV2 
bC  1 I=2,NTftTA 

THLTA(I)=THETA(I-T)4DELT*5.H159/ToC. 

00  7 1 =1  ,NF TiZ 
OC  3 J =1  fNTHT  a 
A =N  4 1 
A X ( A ) = X 

YY(k)=YR(I)*Slk(THETA(J)) 

22U>  = YR(I)*COOUHETa(J)> 

IFLK.Eu* IOC  1 )wk I T E ( NO  , 1 SO) 

IF(A.E<*.10L1)ST«jp 

CON  TINUE 

00  TO  VC G 

CON  T1N0C 

.Ri TE (NO ,2o3  ) 

.81  TL  (NO  ,2  C 4 )L  Abe  L 
.81 TE (NO , It  C ) 

CO  I=1| K 
i Sc  0 = I oFb  + 1 
i F < I* 1 .OT.K  )C0  TO  2 T 
1 F 1 I ♦ 1 • E u • K )E0  TO  c 2 

wRlTE(N0,3CC)XA(l),VY(I  ),ZZ(I)  ,XX(I-H),YV(1*1)  ,ZZ  (I*1),ISE6 
.RiTt(NPt2CC)XX(l),YY(I),ZZ(l),XX(J*1),YY(I4l),ZZ(J*1),lSEfc 

00  TO  ,_C 

• R 1 TL (NO , 3 0 1 )X  X ( 1 ) » Y Y ( I ) * Z Z (I ) ,1lAST,ISEW 
«RlTt(NP,20l)XA(l>,YY(l)fZZ(l)tlLAS7,ISFQ 

00  TO  c 0 

WRlTE(NOt3L2)XX(I)»VV(j)»ZZ(l)»XX(l'*1)1YY(l*1)*ZZ(l41)»ILAST»ISEa 

.RlTt(NP,2C2)XX(l),YY(I),ZZ(I),XX(14l)fVY(l41),ZZ(l41)*ILAST,ISEa 
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tC  CONTINUE 

1 Co  f Oh  HAT  l7A,'A',in,'Y',11X,'Z',13X,'x',11X,'Y',11X,'Z',,m, 

1 / / ) 

15o  F C.>*'  AT  (/ / , 3 > ,'AUhbt  h OF  OFF  tODY  POINTS  EXCEEDS  1000', /,5X, 

1 'FfitCnAh  T£Ff*I  SATED',//) 

F C k K A T (3  ( F10  • 6 ) .1  X , 3 ( FlL«0»1i>X,14) 

,.01  FOKF’AT  (?FlG  «c,l1,<*>X,l4) 

iCi  F0APAT(BF10.t,1X,2Ho.6,11,1AX,l4) 

tC2  I 0 K k A T <1rt1  , SX, 'CALCULATION  OF  THE  OFF  bODY  POINTS  FOR  THE  0CU6LAS 
1 1 Hr  fc  L Olt  ENSIOI.Ai.',/,  OX  , 'POTENTIAL  FLOta  PRCCRAk',//) 
th,  FGr'AT  (S  X.ALA2  ,//  ) 

t L'j  FGrt»ATt4LAc) 

2CC  format c//,e>  .'input  j a t a - location', 1A,//) 

2C  7 FOr^AT  (<>X,'KFTS/-',I**,cX,'NTHETA=',I4,/.6X,'X  = ',F1G»6,2X,'ZSTART- 
1 ,Fl0.c,tX,'tLcYx  = ',FK.t,2X,'TSTART  = 'tFir|.6,2*.'DtLT  = ',M0.6,//) 
:tw  F Ok  *rT  (£  x , 3 (F1l«c,«.X)  f2x,3(F10 .6,2x),1aX,IA) 

- C 1 l Oi\  yAT  (flbiO  iZO  ,1  X,  ll  ,SiA,i4) 

i C c F 0 r ***  r T (2  a t • (Flo»c.t>)»tX,3(FlO»6»2x),1X,l1  .12X.JA) 

SCw  FLr!*AT  (I  j,  1 X ,1  „ oX  , L»F1  C .6  ) 

S TC  F 
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( 

C CALCULATES  TFE  INTERSECTION  POINTS  CF  AN  APPENDAGE  WITH  THE 

c loot  Surface  - apndg3 
c 

S ltd  ROUTINE  AFPINT (X,V,NPTS,  X ST  A R T , X EN  D f LS  E C T , L P T S , X A P P , V A P P , 
1 ZAPP»XINC  ,Y1KC,ZINC) 
b Ih  ENS  ION  X (20u),Y(2t/0> 
dimension  c C EF  (4)  , A (4  ) 

D In  ENS  ION  XINC  (1*.),  Y1NC  (14)  ,ZlNCCU  ) 

DIMENSION  X»FP(1/4,14),YAPK1‘.,14),2APP(l4,l4) 

Din  ENS  IOn  XblRK(5  0,5C;>,Y».C;RK(5u,50>,ZNO«U(50,5i:>,TNETA(50> 

L Tt  ST  = l 
L T = C 
NX=5u 
I T = 5L 

DElT=(  FLOAT  ( IT)  /FLOAT  U T-1)  )*3  .U15W  180. 

THc  TA ( 1 ) = G . 
bO  4 I =2  , IT 

4 THt TA ( I)=Th E TA  (1-1 ) + DEL  1 
DO  5L  L=1*LPTS 

T £i T=C  .0 
Xl  = XAPK1  ,L  ) 

X?=XAPP(LStCT,L) 

Tl-YAPP(lfL) 

Y2=YAPP(L5E  CT,L> 

2l  = 2APM1  ,L  ) 

22-ZAPP(LS£CT,l) 
i 1 L S T = 1 

XSTART=XINC  (t-1) 

CALL  FINEND(>,Y,NPTS,Xl,Yl,Zl,A?,Y2,Z2,XSTAfiT,XEND,YEND,ZEND) 
jO  102  1=1,  NETS 

ICc  I F ( X(  I ).t>E  . XEND  )«G  TU  1 C 1 
1C1  A0t  f=X  (1 -1  ) 

*S T 0P=  X ( I>1  ) 

LTEST=I-5 

LT=I*5 

btu  X = ( xS  TOP-XPeG)  / Fi.0  A T ( NX  - 1 > 

XwORfc(1,1)  = XbEv> 

DO  l 1=2, NX 

3 XnORN (I , 1 )=XwOKX(I -1 , 1 ) 4D ELX 

I/O  5 1=1, NX 
DO  15  .\T  = L T E ST  ,LT 

15  I F (XWOKK  (I  , 1 ).0E. X (NT ) )C0  TO  It 

16  ITt  R = NT~1 

5 CALL  L A6INT  (X,i,PfS  , Y,4  , XtaORK  (I  ,1  ) , tbORK  ( I , 1)  , I IF  R ,1S,C0EF> 

4/C  t 1=1, NX 

DO  t J =1  ,1 T 

AWuRN ( I , J)  = X.OKK(I , 1) 

YWOPK(I,J)=YfcCRK(I,1)*CCS(THLTA(J)) 

( ZWORK ( I , J ) = YWOKM I , 1 )*S IN (THE  TA  (J  ) ) 

91  CONTINUE 

i/O  1u0  I J = 1 , c 
iFOJ.fcU.I  ) IT1=2 
1FCIJ.EU.1 ) 112=20 
1 F(  I J . EQ  .2  ) 1T1  =2u 
1 F(  IJ  . EO  .2  ) I T2  = I T 
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teI3  pack  IS  BBT  QUALITY  FBACTICABLB 

Mom  oopx  tuwishhd  to  ddc~- 


co  1U  l=c,t.X 
kC  1c  J=  111  , lTt 
llM-1 
t.c-  1 
.r  - 1 
.,1=  J-1 

ik  4.  * J “ 1 

iv  5 = J 

AX1=XwCPh(r  1 , N 1 ) 

X At  -XK'Kxtf  t |l<  k ) 

>X^=Xt«Of'A(X  ‘ ,K  j ) 

7 Y 1 -YbChhU  1 ,M  ) 

V Yt  - Yh  OR  X ( l'  k |Kt) 

YY_=Yi»tPx(h  2 ,Ni) 

i Z1  1 ,S1) 

111  = *.Ck^h(hct*kc) 

±2.  - i W 0 >•'  X 

l At  L f L A|  L (>>1  , Y Y 1 , Z* 1 , XX  Z, YVi  , 7iZ , XX  ? , Y Y 5 ,Z  Z 3 , A) 
call  ir.Tsti  oi » y 'i  ,zi,*t  ,y  i ,n  ,a,xt,yt,zt) 

> n t iT -X  x 1 - 7 t S 7 
X ? I LSI =XXt ♦ 7 l S 1 
I 17  1 bT -y Yl ♦ 7 l S T 
»£lli7--YY:-7LSl 
ill lST“?Zl~TtS  1 
* t • l S 7 = Z <.  * ♦ 7 !:  S 1 

.Hx7.t7.xtTiSl.k.f.>l.LT.Al7ES7)tiO  70  cC 
i f l / I 7 _ S 7 . L 7 .I1  )Z1  7 t ST  - _ . 

i F ( >T  . ■ 7 . Y 7 7 ,.S  7 .1  h . Y7  .L  7 . Y^Tfc  S 1)00  TO  tC 

a F l . 1 . i 7 .l  <■  7 t S7  .t  y . i 7 . L 7.  71  7L  S 7 )bO  70  if 

1 Tl  S 7 = - 

XI..  C (L  ) = x7 
VI*.  Ctl  ) - 7 7 
/ 1 ik  t ( L ) = t 7 
l l k 0 . 1 1 N c F 

t F l 17  t S>7  .t  0 .1  ) v.0  7 o Oi 
Ik  ICf.  TINll 

lit  CO  k 7i NOF 

*flI7CS7  >u  •1)7Li7  = Tt.S7'i»;o 
IFIIUL7.LU1)(.0  H r i. 

Oi  C Oik  7 1 h OF. 

!>  L LOi.  7 if.  Of 

h L 7 Oh  l< 
l '.C 


i 
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CONTROLS  THE  C ALL  UL AT  I ON  OF  THt  APP  E ND  A6  E COORDINATES 

i lib  RuUTlNfc  AFPNDt>(X,YtNbODV,Xfc,Hl,NFINtXSTART,XEND,XOFSET,YTlPf 
lYROCT,COKOT,COi<DR,LSiCT,LPTS,XlNC,YINC,ZlNC,XAPPtYAPP,ZAPP,NSECTt 
2 I Ck  Ef  «1C1IP  .TChOOT ) 

DIMENSION  XE  <2b),HC(25> 

DIMENSION  X (t0J),Y (ZOO) 

DIMENSION  XINCCIO.YlNCUOfZlNCtU) 

DIMENSION  XAFP(Ut10,YAPP<U,14),ZAPP(l4,14) 

CALL  F IN  C ku ( Xb  tHl  ,N  F IN  t X S TA  R T *XO  FSET  » YT  I F f Y RO  0 T f C G R D T , C Ok  D R , 

1 LSt  CT  » LPTS  ,X*PF,YAPP,ZAPPfTCREF,TCTlP,TCROOT) 

CAlL  APPINT  (X,  YtiMbODY,  XS  T AR  T , XE  N 0 ,L  SE  C T ,LP  T S ,X  APP  , V APP  , Z A F P ♦ 

1 XIN  CtY  INC* Z INC ) 

NStCT  = LSECT«.i 
DO  27  J=1»LFTS 
xAPPCPSECT  »J)~XInC(J) 

YAFPCMSECT,  J)=YINC  (J) 

L APPCMSECT «J)-ZlNC(J) 

DO  ?0  L = 1 » L F TS 
I F (L .Eu. 1 )OC  To  )u 

DLlY-AdS  (YlNC(L)-YAPP(LStCT  »L)  ) / 2 • 

X I = X AP  P C 1 » L ) 

X2=XAPF(LjECT,L) 

Y1= YAPP ( 1 , L) 

Y2=YAPP(LSt  C T , L ) 

1 V=ZAPr  (1  * L ) 

*2=ZAFP(LSECT,l) 

DO  12  0 I Y = 1 t L 

oO  TO  ( 11 1 1 1 12)  ,IY 
YI».=YAPP(L‘  ECT.L>-DELV 
EO  TO  11i 

Y*i=YAPP<LSECT,l)-2.*0ElY 
CO  NTI MJt 

call  L INECXl  »V  1*21*X2»Y2*Z2  »XFi»YM»Z") 
oO  70(114,112)  , I Y 
XAPPCM  sect-2  ,l)=xm 

YAP  PC M SECT-2 ,L)  = YM 

2 AF  PCPSECT-2  ,L)=ZP 
cO  TO  120 

aAPPCM  SECT-1  ,L  ) =A  P 
YAPPCMSECT-  1 ,l)  = ym 
ZAPPCMSECT-1  ,L 
CON  TINOE 
DO  To  20 

xAPPCMSECT-1  ,L>  = XAFPCLSttT  ,1) 

YAPPCMSECT -1  ,L)=YAPP(LSECT,  1 ) 
lAPP(MSECT-1  ,L)=  ZAPPCL  SECT  ,1  ) 

XAPPCMSECT-*  ,L  ) = XA^PCLStCT  ,1) 

Y AP PCM  SECT-2  ,L)  = YAP PELS  EC T,1) 

ZAPPCMSECT-2  ,L)=  ZAPPCL  SECT  ,1  ) 

COh  T1N0E 

FETLRN 

END 
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1 5o 
5 5 7 


1 I 0 

! 7 7 

i 7c 

j 7 y 


laoMoanf 


^ALI*  PRA.CTICABLB 

tO'CDC  


tSTAcLIShtS  THE  iNPl'T  fCRkAT  Of  THE  CALCULATED  EC  D Y COORDINATES 
fCh  uS£  IN  THE  DuLOLAS  T HK  £ t D IftE N il ONAL  POTE.NTAIL  FLOW  PROGRAM 


SOohullTINE  AAIlu<Xofi>Y,YbODY,ZbCDY 
b Ii*.ENS  ION  XL  LD  Y (IfJO  ,10)  ,YcODV  ClCu, 
LI'ENilOW  X C (1 LOu ) , YL  C1L0C)  ,ZC  C1Q0 
K 0 J M = » 
b C.  1 N - 1 |K 

i f t N.Lc  . 1 ) I STA  lUOJ'.H  1 )=2 
If  ( N . t T . 1 ) liiTrKKOUNM  1):1 
u 0 1 P = 1 , I i 
k Uu  NT  - kO LLTI I •»  1 

aCCKOUNT  ) = X t C D I ( f*  , Xi  ) - X t C D Y ( 1 . 1 ) 
i C ( KoUliT  ) = Y E LC  Y ) 

YC  (KCUM)  = ZL  (0l(i\»fO 

Cui.  TINuE 

iST  AT  (nOuM  ) = 7 

»RiTi:(kQ,3-.-) 

» Si  Tc (£0 , j ^ 1 ) 
i S C.  Jl  — o 

i.!  3 5 5 Ik~1,k0ux1,c 
i!»L>»  = IEt'i'f  I 

ifllMATllK  ).E^.w.«fD.ISTAT(lk-*1). 
ifClSTATdD.M.j.AND.lSTATOk  + l). 
iFCISTAT(lK).E;,,w.«fu.lSTAT(Ik-»1). 
ifdSTiTlH  ).NL.-.hNi>.1STAT  ( 1 k ♦ 1 ) . 
*R1Tl (NO,3t6)XC(ik)  ,YC  (lk),ZC(lk)  * 
00  To  -5a 

*(<ITE(i'.0,3c7)XC(1x),YC(Ik),ZL(iK), 
i zc (i k + i ) , : i lu 

cO  Tu  -51 

of.  1 Tt  (r.O  , Je  • )X  C u k)  , YC  l IK  ) , ZC  ( 1(0  , 
1 I S T A T ( 1 K ♦ 1 ) , 1 A £ <, 

0 C-  T U i 5 5 

(f,C  ,35E)XLClk),YC(Ik),ZC(lY), 
1 «CllK-»1),ISTATClN  + 1),ISEt. 

C C i*  T 1 M c 

if  C lEU.’.Ch.  l u .C  >oG  TO  2 7' 

1 Ll  L-C 

-0  37  5 I f. - 1 t k 0 L f i T ,4. 
iStL=I  LEu-H 

)FClSTAT(l>).Eu.u.AND.lSTAT(lk«1). 
IF(lST»TClk).Nt.L.AfD.ISTAT(lk-»1). 
if  ( IbTxT  (If  ) . E <•  , c .AKo.iSTAT  ( 1 k + 1 ) • 
i M 1ST.  TCIf  ).Nt.u.nM.lSTAT  C 1 K ♦ T ) - 
»«nc(KP,3ct)XCClk),YC(llC),ZC(ik>, 
uC  TO  275 

».  Ri1L<KPf5<7)xC(lK),Y(.(lk)t7t(iK>. 
1 ZL  Uk  ) ,1  5L0 
o 0 T u -75 

* R i Tt  ( k P f 3 f c ) X C ( i k } ,YC  ( 1 k ) , 7 C (IK), 
1 IS  TaT  llkt 1 ) ,1 Stu 
1-0  T u 2 7 5 

*RiTE<NF,30r)XC(itO,YL(IK),7CUK>, 
1 ZC  (Ik  *1  ) , 1 STA  T (1  k+  1 ) , 1 SE  Q 


, k , IZ , I PUNCH  ,NO,NP  ) 
1 C) ,Z60DY( 100.10) 
u ) ,1 STAT (1000) 


EC. 0)60  TO  ?be 
l'.CK-O  TO  357 
UE.OGO  TO  356 
l,  f .0)60  TO  ?3S- 

AC<Ik+1),YC(lk*1),ZC(Ik*1),lSEtt 
1STAT(1K),XC  ( I k + 1 ) .YC(IK+1). 

XCCIk+l)  ,YC(  IK+1).ZC  < I K ♦ 1 ) , 

I STAT(IK),XC  (1  K + 1 ) »YC(IK*1). 


LG.D60  TO  376 
U.^GO  TO  377 
f.  f . C )6  0 TO  3 7 5 
KE.OG0  TO  3 7G 

X C ( I k*  1 ) »YC(  ik*1)»ZC  ( I k-*  1 ) , I S E 0 
ISTAT(lk),XC(lk-»1),YC(lK*1>* 


XC(lk-»1),VC(lK-»1).ZC(IK«1>t 


1STAT(1K).XC(1k«1>vYC(1K«1>* 


.rs 
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sms  HUB  IS  BBST  QUALITY  PRACTICABLE 
*****  00^  HOC  

375  CONTINUE 

35il  FORMAT  ( 2 X,  'FI  NAL  THREE  DIHENSI  ORAL  COORDINATE  OUTPUT',//) 

3 51  FOkMAT  <6X,'X‘,11X,'Y',11X,'Z',6X,'STAT',5X,'X',11X,'Y',11X,'2' 
1oX, 'STAT', 1 2 X, 'St ft',//) 

352  FOKMAT(?X,3(F1C.5,£X),I2,2X,3CF10.5,2X),I2,13X,I4) 

366  F0kMAT(2X,3(F1C»3»2X),4X,3(F1G«5,2X),15X,I4) 

367  F0fcNAT(2X,3(F1o.5.2X),J2,2X,3(F10.5,2X),15X,14) 

366  FOKMAT  t2X,3  CFU.5 ,2X)  ,4X,3(  F 10  ,5,2X),  12,1  3 X,14  ) 

366  FORMAT  <3<F1C,5)*1X*3(F1C,5)  ,15x,14) 

3£  7 F Ok  NAT  <3  < F 1 C.5),1 1 ,?(F 1C. 5) ,1'X,I4) 

3ta  FORMAT  (3(F1C*5),1X,3(F10*5)  ,11  ,14  X,  14) 

3 69  FOkMAT (3 (F1C.5) ,11 ,3(F1C,5) ,11  ,14X,  14) 

3 7u  RETURN 
END 


. * *4+  * *• 

ms  Fiat  IS  MST  QUALITY  PRACTICABLE 
TOM  COPY  PURWISHED  TO  DDC  - 

fALC'LATrS  TIC  BOBy  COORDINATES 

^•nojijTl'JE  DjYAT(X,Y,Nl,TSiXnCDY,YRoDY,7BriDY,NBOOY,!2,NSYM» 

2 T ME  oS  I JN  X ( 2C0  ) , v ( 7 JO  » 

r>  i M e i 5 1 j >1  x n o r>  y ( i j o , i o ) , y b o n y ( i o o » i o ) , ? p n o y ( i o o , i o > 

n t >‘E\'S  1 J'i  TH^T  A I 1 J 1 , COFF ( 1 p 1 , XTFMP  | 1 00  ) , YTE  *P  I 103) 

0 F I = ( * < I **  T S l - X ( l I |/FL0AT(f«'3nr'Y»l  I 

j r ( m -5  Y m . ir  q , i j nr  L T * 3 • M I 59/FL3  A T ( 17-1  1 

! r ( ►'  J Y M , T , 1 ) 01  I.  T s I 9 3 . / F I.  0 A T ( I 2 - 1 >)  • 3 1 B 1 6 9 / 1 e 0 . 

Thrill  1 ) -=0  . 

"i  n 1 1=2,1’ 

! )=THFTA(  I-l  >+'3FtT 

0 F I.  f a f> r L / E>  • 

s"  to  ( n , n , i 3 ) , JSY* 
x r r H r.  ( ( ) _ X ( 1 1 

p o ? I f = ? ( N ° 0 0 Y 

1 ! = I - 1 

0ri.A^3  = O-  Lr*FL0AT(  1 I ) 
ylft"’!  I 1 * X T p M p ( 1-1  1 * 'Y  F L A 0 C 
1 F ( DrLADT. ',F  ,BFi.  ) r,n  TO  27 
L r *'  F>  3 1 

x Tf  I'  3 in  V ) - x ( ,nT  i ) 
n n 7 ) 1 = 7 , ' B 0 n Y 

J '■NsOl:  Y-  I ♦ | 

1 ' = 1 - I 

C F l S ■ ) 3 = Fi  r l F • F L 0 A T ( 1 I > 

v T r >•■  p ( J J ) = X T c "t  p ( J J + 1 I - r F i_  5 !J  B 

I F | 0 F L S J B . 5 F , J E L ) ' 0 T 0 2 S 

,lr*  C SJJ 

JC  tt  0 ] = 1 

l F ( rj$  Y>- . fi  , 3 I if!  T n 3 3 

Y r o s'inoiY-1  F „■  D - j r '•  p I 

rrLKl3=(XTFMP(JFW3)_XTF',B(Lr''0)  )/FL.  FAtIfFIwD^J  ) 
r,  ' 7 rf,  1 = ) , K F <w  0 

J * L F ,v  B + I 

xTr^Bi  jismeiai  j-i  ur-rL'iir' 
g"  To  'O 

Y F V Ds^POBY  - Jr  ivC  I 

rFi  f”  I0=(YTF;*'P(JF>!B|_X(  i i )/float(<f*P) 
xTFKP ( 1 ) = X ( ! ) 
r r 3 7 J = ? , < r ir.  o 
y T F t ' d I 1 ) = XTF>,P  I l-i  ) + 0 c l V I D 
C •"  M 1 M ■ E. 
vTr-'P  i i i =y  i n 
p 0 ? I = 7 , N 3 P u Y 
n r 3 J = 1 , ^ P T 5 

1 F i y r F r P ( I ) , c,  t . X ( j 1 I r' 0 TO  B 
lCyA.?T  = J-  l 

fF,  l-jTAfiT.LF.B)  I G T A S T ■ 1 

C 'LL  L A G j n t | x , M ° T ? , Y ,b,xt£mP(  1 ) ,YTEPT>(  I 1 , I 5 T A F T , I 5 , C C E F ) 
nr  F 1 = 1 , N 3 r 0 Y 
DC  E J = l ,17 
X n n C y ( 1 ,JI*X!EPp<  1 1 

Y Bon  v ( 1 , .))  =YTFw=  ( 1 )«Cr5(TPETA(j)  ) 

7 B nr  y 1 f i J I * Y T F r* p ( 1 1 • 5 1 A ( T H r T A ( j * ) 

K F T U p A' 


»W  WfflfcU  CTM  *UAI<m  PBACTICABLI 


I 

i 


1 

I 

* 


c 

C CALCULATES  TEE  THREE  DIMENSIONAL  COORDINATES  ALONG  THE  PODY 

C 

SUBROUTINE  fc  CD  Y (X  tY , NPT  S , X APP  * TAPP  , Z APP  ,K  , IZ  »M  SECT  »L  PT  S , X B 0 D V , 
1 Y BO  D Y , ZB  OD  Y , XS T AK T f I F * D ,1 Af T ,N S YM ) 

DIMENSION  CCEF(4) 

DIMENSION  XAFP(14,14)fYAPP(l4*14),ZAPP(l4,l4) 

DIMENSION  X E CD Y (100,10) ,VbODy( 100,1 C ) , ZBO D Y < 1 0 L , 10) 

DIMENSION  X1LMF(100),Y1 LMP(10C),THLTA(10) 

DIMENSION  X (20c), Y(2G0) 

N = 1 FWD*LPTS«1AFT 

DEL2=(X(NPTS)-XAPP(MSEC1,LPTS))/FL0AT(1AFT) 

YTLMP (1)=Y ( 1 ) 

CALL  DtLBDY(X(1),XAPP(HSECT,1),XAPP(MSECT,2),XTEMP,IFWD*NSYM) 

DO  2C0  I -2 , 1 END 
DO  1C  1 T — 1 , K FT  S 

1C  IF  (XTEMPC1  ).Gt.X(17))bC  TO  11 

11  1 TE  R=I T-1 

tOO  CALL  LA6INT (>,NPTS , Y , 4 , XTEM P ( 1 ) , YTt PP ( I ) , I TE R , IS , COE F ) 

DO  203  L = 1 , L FT  S 
*TEMP(lFwD+L)=XAPP(MSECT,L) 
c C3  YTEMP(1FWD+L)=YAPPCMSEC1,L) 
laft  = i fwd+lpis-*i 

XTtMP(LAFT )=XAPP(MSECT ,LPTS )*DEL2 
DO  Yt  IT=1,NFTS 

12  IFCXTEMP(LAFT)  .bt.X(IT))GC  TO  13 

13  ITLR=IT-1 

CALL  LAClNT(X,NPTS,Y,4,XTEI*F(LAFT),YTEMP(LAFT),ITfcR,lS,C0EF) 
1C  201  1=2,1AFT 
J=I  F.D+LPTS  M 
XTtMP<  J)=XTF  PP  (J-1)+DEL  2 
DO  14  I T = 1 , N FT  S 

14  iF( XTEPP (J ) .tE  .X(IT  ))GC  TO  15 

15  lTtR=2T-1 

c Cl  CALL  LAG1NT  (X.NPTS,  Y,4  tXTEMr  (J  ) ,YTE  i“P(J)  ,1TER,  IS  ,COE  F) 

DfcLT=(50./FLCAT(IZ-1))*3.14l5V/1cL!. 

THE  TAC1 )=0 . 

DO  1 1=2,12 

1 THtTA(  I)=ThETA(I-1)-»DELT 

DC  2 1=1, 1 FDD 
DC  2 J =1  ,12 
AbODY ( 1, J)  = XTEMF( 1 ) 

YB0PY<1,J)=YTENP(I)*C0S  (THETA  (J)) 

3 ZECDY(1,J)=YTEMP(I)«SIN(THETA(J)) 

i CONTINUE 

X2=  1FV  D*LPT  S -» 1 
DC  4 I =K  2 , K 
DO  5 J =1  ,12 
XB0DY(1,J)-X1EMP(1) 

YB0DY(1,J)= YTEMF(I)*COS  (THE  TA(J)) 

5 2B0DY(1,J)=YTEMP(I)*S1N(THETA(J)) 

**  CONTINUE 

DO  t L =1  ,LP  T S 
1=1  FbD  -»L 

1FCL.NE.DGC  TO  ft 
DC  7 J=1  ,1Z 
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tBIS  Fiat  I MR  practicable 

HUH  OOPT  WBWflHWTO  M>C  " 


1 


aS*0DV<1(J)-JI1EPP(1) 

Y Ho  D V ( 1»J)  = YTEEMJ  ) * 1 0 S (THE  TACJ  )) 

7 ZbO  D V ( 1 »J)  = Y7EnPll  )*SIN(TH£TA(i)> 

bO  TO  t. 

t IOmIIUIE 

YX  = YAFP<MiL  C1,l> 

ZR-ZAPPEEiSE  C 1 * L ) 

K=J>GF\T  (Yk'**?<2M‘J) 

«f»LLL  = ATAN  (Zh/  th) 

i/to  = CV;  . • j,  . 10SV/1f>C«)-ANCLE*i 

liEi.=l»fcv/FLCAKlz-1> 

XbCi)  V (I  ,1  ) = XAFP(KSECT  ,L) 

Ye^OY ( I,1)=YAPF<KStCT  ,L  ) 
iPOCY  < I .1  ) = ZAPP(i*iStCY  , L ) 

Abb  D Y (i«IZ  > = X A F P l K b F t T , L ) 

Y Pw>  D Y (i*lZ)  = ZAHP(HS(CT,L) 
ibl)l»Y  (J  ♦ I Z ) = YAbF(MirCT  »L) 

L2- 12-1 

J Ft  TA ( 1 ) =AN£Lt 

I/O  ■>  J=?,LZ 

lHtTA(J)=ThfclAU-1)-U>tL 
XowDY<i,a)  = X1EKrU) 

Y t o t)  Y (1  *J)=-P»CGSlTrtt1A  (J)  ) 
i?tDY(l,J)  = H*$l.\(TH£TA<J>> 

S-  lOr<T]\t,£ 

■j  IOM1ME 

ELI  GEN 

L b 


snsstfs.: 

fKMOOft] 


,ln  KRACX1CABU 

>£Jbc 


l 

C CALCULATES  TEE  SPAC1N6  OF  THE  POINTS  ALONG  THE  BODY  SURFACE  IN  THE 

C X 0 1RECT10N 

C 

SUd ROUTINE  D tL  bD V (X  *XA  PP1  , X A P P 2 ,XTE  PP  * I F WO  ,N  S YM  ) 

DIMENSION  XTEMF(IOC) 
iF(lFWD.LE.6)GC  TO  40 
DEL =<X APPl-X)/ FLOAT <1  FDD) 

DELFlN=XAPPt-XAPPl 
00  TO  <21 ,21 t2i>, NSYM 

21  XTt  MP ( 1 ) =X 

DO  1 I=2*IFRC 
11=1-1 

DElADD=DELF  IMFLuAT  (II  ) 

ATEMP<l)=XTEkP(l-1)*DELADD 

1 iF(DELADD.6E.DEL)G0  TO  2 

2 LfwD-I 

22  XTEMP<IF«D>=XAPP1-DtLFlN 
DO  4 1 =2  *1 F Db 
JJ=IFHD-1*1 

DEL  SUB=DtL  F Ik*  FLOAT  <1 ) 

XTtMP(JJ )=XTEMPlJJ*1)-DELSUB 

4 1 F ( Dfc L SUb. 6 E .DEL) 60  TO  5 

5 JFhD  = JJ 
JF*DI=I 

1F(NSYM.EU.3)6C  TO  10 
KF*D=1 Fub-LFxD-JFHDI 

DELPID  = (XTEkK  jFhDI-XTEPP  <L  F W D ) ) / F L CAT  (X  F W D*  1 ) 

DO  t I=1,KFWD 
J =L  FN  D ♦ I 

t XT2HK  J)=XTEFP(J-1)+DcLMl  D 

RETURN 

1C  xFb D=I FND-J FhD  1 

DELPlD=<XTEkF(JFfcD)-X)/FLOAT(KFWD> 

XTE  MP ( 1 ) =X 
DO  12  1=2,AF*D 

12  XTEMp (1)=XT  E PP  (1-1 ) *D  t L HI D 
RETURN 

-0  DEL=(XAPP1 -X )/ FLOAT <1 FW D-1) 

X TE  MP  < 1 ) =X 
DO  A 1 1 = 2*1  FDD 

41  XTtPP<l)=XTEFP(I-1)*DEL 

RET  URN 
L NO 
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1H23PAQK  IS  BEST  QUALITY  PRACTICABLE 
JROM  OQPY  flBaisrifiD  TO  OCC  „ 


1 1 


-i59  ki<llL(KO,3^)l(C(n),VC(lK),K(lK)(ISTAI(lK),XC(IKO),YC(IKM)f 
1 iC(IK+1),IST/»UU+1),lS[w 
j55  CONTINUE 

IF(IFUnCm.EU.0i6U  TC  3 7 C 
i S e.  u = c 

i,C  375  I k-  1 ,XOCNT,c 
ISlGM  SE«*1 

IF  ( 1ST  AT  (IN  ) .E  «..u.Aki>.l  ST  AT  (lk«1)  .bfl.OGO  TO  376 
1 F ( 1ST  AT  II K ) .\t . J.AND  .1 STAT  (I  X ♦ 1 ) .E  3.0)60  TO  3 77 
1ST *T (I N ) .E *.j.AN0 .1  ST  AT <IK*1) .n E .0)60  TO  37? 

1 F ( 1ST  AT  (IN  ) ,M  .j.mM,  .1  ST  AT  < I K •»  1 ) .N  t .060  TO  379 

.PITt<NP,jt6)XC(ilO,YC(lK.),2C(IIO,XC(IlC»1),YC(iK-H),ZCCIK'*1>,ISE* 

00  TO  Z^S 

3 77  KKl  Tl  (NP  ,'ii>  7 )Xl  (IN  ),YC(IlO,ZC(lK),lSTAT(IK),XC(IK*1)vVC<IK-»1>, 

i :c  (1k+i  >,  i stci 

»0  TO  375 

-7o  iti  Tc  ( N P , 3 <.  £ )Xt(iK)«YC(IX  ),ZC(ik),XC(IK+1),YC(lK+1),ZC(IK+1)t 
1 IS  TAT  (Ilk+  1 ) ,1  itu 

00  TO  373 

3 7 v wk.tTE(NP,3E9)Xtax),YC(lk),ZC(lK),lSTAT<IK),XC(IK+1),YC<IK  + 1), 

1 Zt  UK  *1  ) , 1 STAY  (i  «♦  1 ) , 1 St  U 

2 7 3 CONTINUE 

j5u  F 0 *\  F a T ( 2X,  'FINAL  T m Ft  E D 1 k L N S I u N AL  COORDINATE  OUTPUT",//) 

351  F 0 k F A T 16  A, " X " , 1 1 X » " Y " , 1 IX ,"Z",tX,"STAT",5X ,"X",11X*"V",11X  * " Z " » 

1 i X f STmT  , 1 £ X , s l 0 , / / ) 

^52  F 0 r F a T (2x,^<f1c.3»uX),l£t£X,3(FlG.5»2X),l2,13x*14) 

3 1 c FGkFATt2x,3<no.3,2X),4X,?(F1C.5,2X),15X,I4) 

367  f0k?AT(ZA,3OlLO,cv),i2t2X,3(f1u.5,2X),15X,14) 

3iu  FOG  FAT  (fA  , 2 (Hc.j,2>),4X,3(F1C.5,2X),I2,13x,14) 

3F6  F Or  F A T ( 1 (F1C.5)*1X,71F1L.5)  *15*, 1**) 
if  7 FCkFaT  < 3 ( F 1 C . 5 ) , i 1 ,2(F1C,5)  ,15  A, 14) 

;id  FOrF AT  <3 (F 1 C .5 ) ,1 X , 3 (Fl 0.5) ,11 ,14X,  14) 

:,lv  FO  FF-AT  (3  (FI  C .5)  ,i1  ,3  ( F 1 C.  5)  ,11  ,UX,  14) 
i 7o  htTL'RN 

t N L 


aBiamoKia  awaQumm  practicabu 


c 

( LSI AdL  1SHE  S THc  INPUT  FORMAT  OF  THE  BODY  AND  APPENDAGE  COORDINATES 

C FOR  USE  IN  THE  DOUGLASS  THREE  DIMENSIONAL  POTENTIAL  PROGRAM 

C 

SUBROUTINE  DCU6C  (X bOO Y ,Y bO D 1 , Z BO D T ,IC , 1 Z , X APP , » AP P , Z APP  ,MS E CT  . 
iLPT  S , I PUNCH  , NO  ,NP  ) 

DIMENS  ION  XAFP  (U,14)  ,YAPP(14,1A),Z  APP  (14.14) 

DIMENSION  A6CDt(100,10),YeODY(100,1C),ZeODY(10C,10) 

DIMENSION  XC<1Sao)*VC<15C0)  ,Z C ( 1 5 OG ) , I ST AT <1 50C > 

KOU  NT  = G 
DO  1 N=1,K 

1 F(N.ED.I)  1 STATUOUNTM  )=2 
IF (N.CT.1 ) ISTnT(KOUNT41)=1 
CO  1 H=1  ,12 
COUNT =KOUNT  -*1 

XC  ((COUNT  )=XECDY(N,K)-XbODY(1,1  ) 
t.  C ( (COUNT  )=YUCD Y(N,M) 

YC((COUNT)  = 2E-CdY(N,M) 

1 CONTINUE  I 

i/O  2 N=1  ,LP1S 

IF  (N.EG.1  ) 1 STAT  OcOUNT+1  )=2 
1 F (N.  61.1)1  STA  TUOUNT  + 1 )=1 
i.C  2 "=1  »MS  E CT 
f,M  = MSE  CT-M*  1 
kOuNT  = NOUST  *1 

*C(KOUNT)  = X AFP(MntN)-XbCDY(  1,1) 

/ C ( KOUNT ) = 2 A FP  (MM , N ) 

Y C ( (COUNT  ) = Y A FP  (KM  , N ) 
c (ONT1NUE 

DO  3 N =1  , LP T S 

IF  (N.EQ.1  ) I ST  AT  ((COUNT-*  1)=2 
1F(n.GT.1)IS1A1(|COUNT«1)=1 
DU  3 M =1 ,M$  t CT 
NOUNT=KOUNT  •*  1 

XC(KGUNT)=XAFP(M,N)-XdODY(1  ,1) 

2C(ICOUNT)  = VAFP  (M,N) 

YC(ICOUNT)  = 2APP(M,N) 

3 CONTINUE 

I ST  AT  ( (COUNT  ) =3 
• RITE  (NO, 350 
bRl TE (NO  ,35  1 ) 

1 St G-C 

DO  355  IK  = 1 ,KOONl  ,2 
!StU  = I SE  U+  1 

1F( JSTAT(lK).EC.u.AND.ISTAT(IK*1).EG.0)GO  TO  356 
1 F C ISTAT  (IX  ) .Nt.U.ANO.I  STAT  (I(C-*1)»£'3»0)60  TO  357 
I F( 1STXT (I  X ) .ED.Q.mND  .1  ST  AT (IX41)  ,NE.C)60  TO  356 
IF(ISTAT(IN).NE.U.  AND.  1 ST  AT  ( I K + 1 ) . E . 0)60  TO  359 
3 56  *RITE(NO,3o6)Xl(IIC),YC(llc),2C(IIC),XC(llC-»1),VC(lK-»1),ZC(IK+1),ISEt 

GO  TO  355 

357  •kiTt(NO,3o7)XC(itC),rC(IIC),ZC(IK),ISTATClK),XC(IKA1>,YC(IK4l), 

1 ZL(I(C+1),ISFG 
oO  TO  355 

3 58  *Ri  Tt  (N0,36F  )XC  (IX)  ,YC  (l(c),ZC  (IK)  ,XC(IK*1  ) ,YC<  IK-M  ),ZC  ( I (C  ♦ 1 ) , 

1 IS  TAT  (IK«1  ) ,1 SEG 

GO  TO  355 


SHIS  PKtt  IS  HUM  QlttMTY  PRACTICABLE 
ISOM  OOPTf  fUfWISHED  K>  DDC 


C CAL  lULATtS  7FE  FIN  ClGkOINATtS  FROM  THE  GIVEN  INPUT  COORDINATES 

C hHii  THt  STARTING  POINT,  CORO  VALUES,  AND  OFFSET 

C 

SUL  ROUTINE  FiKLnti(XP,HC  ,1PTS,XSTART  ,XOFSET  ,YTIF,YROOT,CORDT,CORDR, 
IlSlCI  » L P T S , XxPP,YAPP,aAFP,TCREF,TCTIP,TCROCT) 
uIi'iENSiOi,  C C t F ( <i ) 

tirUKSiCN  ).  A|PCn,14),YoFr(lA,14),2TPPC1A,14) 

DIMENSION  Xf(2S),HC<c5>) 
i £l=1 . 1/ FLO  A T ( LP1 S-2) 
u 0 1 I -1  ,L  S f CT 
A APT  (1,1 >-L  . 

A A p P ( I ,Z)-UEL/t. 

A AF  P(I  ,3  ) = 0 F.L 
i/O  1 J-A,LPTS 

I XAPP(I  ,J  )-AAFP(I  ,J-1>*DtL 
i/O  v I =1  ,2 

GO  TO  (3  ,4  ) , 1 
0 L - 1 

o(j  5 J -1  ,LP  1 S 

UlL  LAG  iNT  (*P  , I r’  T S , h C ,A,XAFP(L,J  ),ZAPP(L,J)  ,1  ,IS,COEF) 

aAPP(L,J)=XCFSuTAXSThRT+*APP(L,J)*CCkDT 

lAt-  r (L  ,J  > = Y T IP 

b :APP(L,J)=iAIP(L,J)AlORoT*(TCTiP/TCFEF) 

a tc  ' 
c =l  S L C T 
u 0 P j-1,LPTi 

LAgL  LAC  i N 1 OF  ,IF  TS  ,hC  XAi  P(L,J  ) , 2APPCL  , J)  ,1  .1  S ,COEF) 
a Ar  P ( L ,J)  = aSTAkToAPP(L  ,J  OCORDF 
Y AfP(L  ,J  > = Y F iO 1 

t 2 Ar  P(L  ,J  ) = *'  a FP  <L  , J > *COhOK*(  TCROOT/TCkEF) 

. l On TINgE 

1/ElFiN-<yt  i f - y r gg  t ) / float  (lsect-i) 

i,?L  CT  = LSlCT  - 1 
.0  1w  I-2,KSt.CT 
v 0 10  J = 1 , L * 1 S 

1C  Y 4 P F d »J)-YA*P(I-1tJ)-0fcLFIF 

l 0 11  1 = 2 * P SECT 

i V 11  j=1,lFTS 
Xl  = AhPF ( 1 , J ) 

Ai  = XAPP(LSr.  C1,J) 

Yl=YAFF(1,j) 

Y t = YAPKLS  L C T , A) 

J-'  ZAP  F ( 1 , J ) 

,2=ZAPP(LSF  C1,J) 

II  UlL  L IN t < A 1 ,Y 1 ,l 1 , XL, Y 2, 22 .X  A rP (I , J) , YAPP  (I ,J  ),ZAPP U ,J)  > 

RET  D R N 

Lot/ 
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9ns  PUB  IS  NBSJSWtf™  PRACTICABLE 
mOM  OOPY  IUHBISKED  TO  DDC 


C 

C CALCULATES  THE  APP  POX  I FATE  INTERSECTION  PCINT  OF  THE  APPENDAGE 

C TRAILING  EDGE 

C 

SUb ROUTINE  FIN END (X ,Y ,N60DY ,X1 , Y1 , Z 1 , X 2 * Y 2 ,2 2 , XST ART , X EN D , YE ND  , 
1 ZEND) 

DIMENSION  C C EF  (4) 

DIMENSION  X <200,  YC200) 

1ST ART=0 
DO  1 1=1 tNBCDY 
1F(ISTART.G1.0)GU  TO  1 
IF(X(I ).GE .XST ART )1 STAR T= 1-1 

1 CONTINUE 
D EL  = • 0 05 
TES  T=.G5 
1TER=0 
1M=  Y2 

2 I TE  R= I TE  R*  1 
Y N = YM- DEL 

CALL  LINE(X1,V1,Z1,X2,Y2»Z2fXM,VM,ZP> 

CALL  LAG1NT  (X,f»bODYtY,4*X*,YMl  ,ISTART,IS,CCEF) 

XTtST=ABS(YP-Yh1) 
iF(XTEST.Lt.TEST)60  TO  5 
IF(ITER.6T.5C0u)iiO  TO  5 
GO  TO  2 
5 X EN  D = X P 

YEND=YM 
ZEN  D = ZM 
RET  URN 
END 


! 


SBIS  PAG®  IS  BBSS  QUALITY  PJU.CIICABI£ 
TBM  COPY  rURBISHH)  TO  DOC 


l 

t L A»_  CliL  AT  t S TFE  I *T  t F S t C.  T I bN  POINT  OF  A THRfcE  DIMENSIONAL  LINfc 

C ANu  A JHKfct  PINEhSlCNAL  PLANE 

C 

c 

C tibATlCN  Of  LIhc  (X-X1 >/ L=  ( V - Y 1 > /M  = ( Z -Z  1 )/N 

c l=iosl  *=cos*  n=ccsn 

t 

• l tUciATICN  OF  FLANk  A » X ♦ c*  T+  C * 1 ♦ D = t 

L 

itibFOUTl.it  lNTSEL(Xl,Y1t?1,XitVZ,Z2'C0tF(Xl,YT,ZT> 

. DIMENSION  CCtFU) 

LD-S£iRT( ( X 2 - Xl )**2+(Y2-Yl )**Z+(Zc-Zl)**2) 

tOSL=(x2-Xl  ) y D D 

LOS*=( Y ? - Y 1 ) / 0 0 

L0oN=(t2-Zl  >/Du 

*=CCEF  Cl  ) 

c=cCtF  (?) 

t = C 0 k F (3) 

i-=CCEF  (A) 

I F ( X2  . t3  .A  T KO  1 ^ 1 

aT=<(i.  *(COsi*/CCS».)  + C*((CSN/tOSL))*Xl-D-P*Y1-C*Zl>/(A4B*(COSN/COSL> 
1 ■*  C * (CC'SN/Cc5l)  ) 

Y T = Y 1 4 ICOSN  /LGbL)*(XT-X  1) 

4.T-Z14  (CJSN/t0iL>4(XT-Al) 

NET  UnN 

1 C On  TIN  OF 

Y T-(<A«(C0LL/Ct5h>4C»(CGSN /COS  O)*Y1-A*X1-C*Z1-D)/(A«(t0SL  ZCCS  H)*B 
1 ♦ C * (C3SN/C0S>)  ) 

Al=xl4(CtSL/l0a*>*(YT-Yl> 

X T = 1 1 ♦ (CoSN  / LO  SM ) * (YT-V  1) 
k t T Ij  K N 
L W 
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SHIS  BASE  IS  BEST  QUALITY  PBACTICABLB 
IStOtf  OOIPY  FUBKISHED  TO  


SUBROUTINE  LAC1NT  < VI N T ,N 1 , Vd E P ,NPT S , V IN , P S , J , I S , COE F > 

C 

c 

C OBTAINED  FRCP  D.M.  NELSON  NOSC 

C 

C THIS  SUBROUTINE  PERFORMS  LAB  RANGE  INTERPOLATION  FOR  ANY  DECREE 
C LESS  THAN  A 

C VINT  = TABLE  OF  VALUES  CF  INDEPENDENT  VARIABLE 

C NI  - NUMBER  OF  ENTRIES  IN  VINT  AND  VDEP  TABLES 

C VDEP  = TABLE  OF  VALUES  CF  DEPENDENT  VARIABLE 

( NPTS  = NUMBER  OF  POINTS  FIT  EV  POLYNOMIAL  (DECREE  PLUS  ONE) 

C WIN  = VALUE  OF  INDEPENDENT  VARIABLE  TO  BE  MATCHED 

C PS  - MATCH  VALUE  CF  DEPENDENT  VARIABLE 

C J - INDEX  AT  HHICH  TC  START  StARCH 

C IS  = SUBSCRIPT  OF  FIRST  POINT  OF  THOSE  FIT  BY  POLYNOMIAL 

C COEF  - ARRAY  Of  LAC  RANGE  COEFFICIENTS  USED  IN  EVALUATION  OF 

C DEPENDENT  VARIABLE  ; PS(VIN)  = COE  F (1  ) * V D EP  ( I S ) ♦ *** 

C 4 COE  F(NPTS) *VDE  F( IS+NPTS-1  ) 

C * ; 

fc  FORMAT  (/2 X ,58HLm6I NT  EXTRAPOLATEO  MATCH  VALUE  FOR  INDEPENDENT  VAR 
11ABLE  -*£12  mt/ ) 

7 FORMAT  </2X  .IOSHlAoRANCE  INTERPOLATION  (SUBROUTINE  LAG1NT)  ATTEMPT 
1 1 D FOR  TOO  IaRoE  A Dt6REE  - PROGRAM  TERMINATED  * NPTS  = ,I3) 

DIMENSION  XT  U),YTU>,VINT(N1),VDEP(N1), COEF  (NPTS) 

IF  (NPTS.GT.O  60  TO  60 
IF  (J.LE  .1  ) J = 2 
IF  (J . CT  .NI  ) J = NI 
IF  ( V I NT ( 1 ) ~ ViN)  25,51,97 
25  IF  (VINT  (NI)-  VIN)  96*51,51 

51  IF  (VINT  (J  )-VIN)  52  ,56,56 

52  J = J+1 
CO  TO  51 

56  IF  (VINT  (J-D-VIn)  2o  ,2  2,57 

57  J = J-1 
bO  TO  51 

it.  IF  (J  - NPTS72  ) 96  ,96,29 

29  IF  (J  ♦ NPTS/2  - NI  ) 3C, 30 ,95 

30  IS  = J - NPTS/c 
33  IE  = IS4NP1S-1 

L = 0 

DC  31  K=IS,IE 
L = L4l 

A T(  L)  = VINT  (K  ) 

YT(L)  = VDEF(K) 

31  CONTINUE 
PS  = 0.0 

DO  AO  M=1,NPTS 
PF  = 1.0 
DO  29  N=1,NF1S 
IF  (N-M)  32  ,29,32 

32  PF  = ((VIN  -XT(N))/(XT(M)-  XT(n>))«  PF 
39  CONTINUE 

COfcFiM)  = PF 
PS  = PS  4 Y T (M ) * PF 
AO  CONTINUE 
CO  TO  9 A 
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9t  wRITt  (6,6)  VIN 
VS  IS  = N1  - NFTS  ♦ 1 
00  TO  33 

•i  1 wRl  Tt  (6  ,6  ) VI  is 
H iS  = 1 
CO  TO  -3 
*)<«  KtTl'RK. 

60  WRITE  < 6 , / ) NPTS 
RfcTl'K* 


* - 

Ibis  pa ge  is  best  quality  practicable 

FROM  COPY  FURNISHED  TO  DDC  -- 


.4W 


HBIS  PAGE  IS  BIST  QUALITY  PRACTICABLE 
JBOi  OOPI  FURKISHESD  TOLDC 


C 

C CALCULATES  TEE  PARAMETERS  OF  A TMREEE  01  RE  NS  I ON AL  LINE 

C 

SUBROUTINE  L IN E (X 1 , Y 1 , 2 1 , X2  * V2 ,Z2 , X P * VN* ZN ) 
D=SRRT(CX2-X1)**2ECY2-Y1)**2+(Z2-Z1  >**2> 

COS  L* ( X2-X 1 ) /O 
COSN=< Y2-Y1 >/D 
C0SN-CZ2-Z 1 I/O 
XH=  X1  + (C0SLSC0SN)*(YM-V 1) 

ZN*Z1*  (C0SNZC0SNWYH-Y1) 

RETURN 


Ibis  page  is  best  quality  practicable 

FROM  COPY  FURNISHED  TO  DDC 


COMPUTES  THE  COEFFICIENTS  FOR  THE  EQUATION  OF  A THREE 
DIMENSIONAL  FLANt  l»  1 W E h THREE  POINTS  ON  THE  PLANE 


LQOATItN  OF  FLANt  C OE  F (1  ) * X*  C 0 E F ( 2 ) * Y *C  0 E F C 3 ) ‘24  COE  F < 4 ) =0 

sUt  ROUTINE  FLAhElXl  ,Y1  ,Z1  fX2,Y2  ,Z2,*3t  V3.Z3.COEF) 

DIMENSION  CCLFU) 

DtT  = Xl  *(Y2*Z2-Y3*Z2  )-xZ*(Y1*Z3-Y3*2  1)-»X3*(Y1*Z  2-Y  2*Z  1) 
Al=-<(Y2*Z3-Y3*Zc>-(Yl*z3-Y3*2l>+<Yl*Z2-Y2*Zl>) 

A.  2-(x2*Z3-X3*Z<.)-(Xl*Z.3-X3*Zl)-*(Xl*Z2-X2*Zl) 
a3  = -((x2*Y3-x3*Y2)-(x1*Y3-x3*Y1)'*(X1*Y2-X2*Y1)) 

COt  F(1 ) = A1 / DLT 
cot  F t2 )=A2 / D ET 
tOt F (3 ) = A 3 / DET 
c OE  F ( *»  ) = 1 • 

RET  CRN 

tf.D 
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BX9  NWttJS  BUST  QUALITY  PRAC1ICABL? 
jbom  <xipy  Vurjlished  ro  ucc  — -- 


c 

C PLOTS  THE  CALCULATED  THREE  DlMtNSlONAL  COORDINATES  FOR  AN 

C AXISYMMETRIC  BODY  BITHOUT  APPEUDAbES 

C 

SUBROUTINE  FLOTACXbODY  tYBODYv2bODY  , N P T S » I Z ) 

DIHENS  ION  XeCDYU00,10),Yfe0DY(100*1C),ZB0D\C10C,10) 

DIMENSION  X PLO T O OU ), YP LO T ( 100 ),Z PLOT C 100) 

1ST  ART  =0 

YU--ABS  ( YBOC  V 1 1 ,1  ) - Vb  3 D Y C 2 , 1 > ) 
i.  0=  YO 

AHA  X = XbO  D Y CRFTS.1 > 

YNA X = 0 . 

DO  70  1=1 fNFTS 

7 C iFCYBODYCI ,1  ). bT. YM A X ) Y HA X= YbO 0 Y ( I ,1) 

*MAX  = VMX 
1 PL  0T  = 1 
M = J CO  » 

D=S0C. 

c=soo. 

*C  = XbODY  (1  , 1 1-AbS (XBOD  Y (2  ,1 )-X60D Y ( 1,  1 ) ) 

CALL  BGNPL  (1FLCT) 

CALL  PAGE(15.«11.) 

CALL  TITL3D  ('  4 1 00, 1 0 . , 10  . ) 

CAlL  FRAME 

CALL  AXES3D<*X»',1o0,'Y4',1C0.'Z>',1C0t0.,b.,0.) 

CALL  VOAbS  (A  ,B  t C ) 

CA»_L  feRAF3fc(XOt'SCALE'tXMAX,YCt'SCALE'’,Y«Ax,ZOt'SCALE',2MAX) 
DC  60  I=1,\P1S 
DO  61  J = 1 f I Z 
XPLOT (J)=XBCDY  ( I , J ) 

YPLOT ( J)=YuCDV  (1 , J ) 
tl  aPlOT ( J)=Zt CDYlI, J) 

6l  C Ai.  L CURV3C  (xPLCI»2PL0T,YPLCTtIZt0) 

DO  63  J = 1,  1Z 
DO  tb  1=1,NFTS 
XPLOT (l)=XbCDYU.J> 

YPlOT C I)=YBCCY  (I , J) 

66  ZPLOT ( i)=ZoCDY ( I » J ) 

0 5 C Ai.  L C JR«/3D  OPLOI  ,Z  PLOT  , YPLOT  ,NPTS  , D 

CALL  ENDPL(IFLLT) 

110  CONTINUE 
Ret  URN 
t NR 
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c 

C PLOTS  THL  CAICLLmTED  THREE  DIKtNSlONAL  COORDINATES  FOR  A BODY 

L »ITH  APPENDACES 

C 

bUr,mUTIi«E  PLOTCCXbCOY  ,YbCDY,7oCDY  ,NPTS,1Z  ,XAPF,YAPP  ,2APP,NSECT, 
1 LPT  S) 

l IwENS  ION  x»FP(1<.,14),YAPP(n(l4>,2APP(14,14) 

DIMENSION  xECDYn0u.10),Yb0DY(V30,1C)t2B0DY(10C,1D) 
l It.  INS  ION  XFLOT  C 1 Go  ) , YPLOT(  1 tiQ)  *2  PLOT  ( 10L  ) 

1 ST  ART  =0 

YO--Ac l < YAP F 11  ,1)-YAPP(c,1>> 

Z 0 - Y C 

xw«x-xl01/t  ( a fts  ,1 ) 

Yr-A  X = Y APP  ( 1 ,1) 

Z!*i*  a * Y M X 
m =i  CU . 

L “J  Cu  • 

C=JnC. 

1 PL  GT  = 1 

xO-XbODY (1 , 1)~ApS (XbODY  (2,1 )-XoOCY(1,1)) 

Ha  BC-NFL  (lPLuT) 
call  f AGt (1 5 .,11. ) 

CALL  T 1 T Li  C (*  i 1LC, 1C  ., 1u  .) 

CAlL  FrAKE 

UlL  AXESjD(*Xi',1uC,'YS',10C,'ZJ',1uC,0.,L.,0.) 

CAlL  VCiAt>S(A,b,C) 

l A..  L tRAF3L<XO,'SCALt',XF,AX,YG,'SCALE',YnAX,ZO,'SCALE',ZPAX) 
c 0 ;l  1*1 ,P  StC  T 
CO  il  J=1,LPTS 
xFLOT ( J)-X  A F F ( 1 ,J  ) 

YPLCT  <J)=YAFF  ( i , J ) 
i 1 lPlCT  (J)=ZAFF(1,j) 

CAlL  CuRVil  OPLOf  ,YPL0T  ,2PLCT,lPTS,C) 

3[  CALL  CJRV3D  (XFLOT, ZPLOT, YPLOT  ,LPTS  ,C) 

t-U  iu  J=1,LFTS 
CO  M 1 = 1,  A SECT 
XFLOT  ( 1 ) =X  A P P ( 1 ,J  ) 

YPLOT  ( 1 ) =Y  A F F ( i ,J  ) 

Cl  ZPLCT  ( 1)=ZAFF(  . , J ) 

LAlL  C u R V i C (aPlOI  , YPLOT  ,ZFL0T  ,i*  SECT  » G ) 

CALL  CcRViC  OPlOT  ,Z  Plot  , YPLOT, (1SECT  ,C) 
u C Cu  i*1,nfts 

00  ( 1 J=  1 , I z 

XPLCT ( J)=Xt ( DV  (1, J) 

YPLCT ( J ) = Yt  CC V (I  , J ) 
t 1 JcPT  <J)=Zk  CtYCi.J) 

6C  CAlL  CCIPV3C  ( aPL^/1  , Z PL  0 T ,YFLGT,IZ,0) 
lC  cb  J=1,1Z 
LC  f t 1=1 ,NPTS 
>PlCT(  1 ) = X l IDYU.J) 

YPLCT  ( 1 ) =YF  C C Y (I  , J ) 
it  .Plot  ( I)=2kC0Y  u,j> 

c 5 LAlL  ccrvjl  oplgt, zplot  ,vflct,npts  ,C> 

LAlL  ENDFL ( IFLlT) 

L A,.  L t ONEPL 
RLT  URN 

1 ND 


83 


